3Y,  CA-.  :.A  96843-6004 


NAVAL  POSTGRADUATE  SCHOOL 

Monterey,  California 


^  (o(-5'/a 


COMPUTER  ALGORITHMS  FOR  MEASUREMENT       | 

CONTROL  AND  SIGNAL  PROCESSING 

OF 

TRANSIENT  SCATTERING  SIGNATURES 

by 

Soonpuen  Sompaee 
September  1988 

Th 

esis 

Advisor:            Michael  A.  Morgan 

Approved  for  public  release;  distribution  is  unlimited 


T242369 


iclassified 


jnty  classification  of  this 

page 

REPORT  DOCUMENTATION   PAGE 

Report  Security  Classification  Unclassified 

lb  Restrictive  Markings 

Security  Classification  Authority 

3  Distribution  Availability  of  Report 

Approved  for  public  release;  distribution  is  unlimited. 

Declassification  Downgrading  Schedule 

Performing  Organization  Report  Number(s) 

5  Monitoring  Organization  Report  Number(s) 

Name  of  Performing  Organization 
aval  Postgraduate  School 

6b  Office  Symbol 
(if  applicable)  32 

7a  Name  of  .Monitoring  Organization 
Naval  Postgraduate  School 

Address  (city,  state,  and  ZIP  code) 
onterey,  CA  93943-5000 

7b  Address  (dry,  state,  and  ZIP  code) 
Monterey,  CA  93943-5000 

Name  of  Funding  Sponsoring  Organization 

8b  Office  Symbol 
(if  applicable) 

9  Procurement  Instrument  Identification  .Number 

Address  (city,  state,  and  ZIP  code) 

10  Source  of  Funding  Numbers 

Program  Element  No    Project  .No     Task  .No     Work  Lnii  .'Accession  .No 

Title  (include  security  classification) 

:OMPUTER  ALGORITHMS  FOR  MEASUREMENT 

:ONTROL  AND  SIGNAL  PROCESSING 

)F  TRANSIENT  SCATTERING  SIGNATURES 

Personal  Auihor(s)  Soonpuen  Somapee 

a  Type  of  Report 
aster's  Thesis 

13b  Time  Covered 
From                        To 

14  Date  of  Report  (year,  month,  day) 

September  1988 

1 5  Page  Count 
91 

Supplementary  Notation  The  views  expressed  in  this  thesis  are  those  of  the  author  and  do  not  reflect  the  official  polic>  or  po- 
ion  of  the  Department  of  Defense  or  the  U.S.  Government. 

Cosati  Codes 

18  Subject  Terms  (continue  on  reverse  if  necessary  and  identify  by  block  number) 

Radar  Cross  Section,  Transient,  Scattering,  Signature 

;id 

Group 

Subgroup 

Abstract  (continue  on  reverse  if  necessary  and  identify  by  block  number) 

This  thesis  describes  the  development  of  computer  algorithms  for  experimental  measurement  control  and  subsequent 
mal  processing  of  transient  signatures  to  synthesize  scattering  impulse  responses  of  scale  model  targets.  The  theories  behind 
msient  scattering  are  considered  in  order  to  construct  the  algorithms.    Up-to-date  hardware  and  software  technology  are 
lectively  implemented  to  optimize  the  resultant  signal  to  noise  ratio.  The  detailed  explanation  of  the  hardwiu-e  and  software 
)eration  is  provided.    A  noise  model  for  the  system  is  also  discussed.    Measurement  and  programming  validations  are  de- 
ribed,  where  comparisons  are  made  with  numerical  computations  for  transient  scattering  by  selected  canonical  targets, 
esults  are  documented  for  futher  research  in  the  transient  scattering  problem. 

Distribution  Availability  of  Abstract 

unclassified  unlimited          D  same  as  report          D  DTIC  users 

21  Abstract  Security  Classification 
Unclassified 

a  .Name  of  Responsible  Individual 
Michael  A.  .Morgan 

22b  Telephone  (include  Area  code) 
(408)  646-2677 

22c  Office  Symbol 
62.\Iw 

)  FORM  1473,84  .MAR 


83  APR  edition  may  be  used  until  exhausted 
All  other  editions  are  obsolete 


security  classification  of  this  page 
Unclassified 


Approved  for  public  release;  distribution  is  unlimited. 

Computer  Algorithms  for  Measurement 

Control  and  Signal  Processing 

of  Transient  Scattering  Signatures 

by 

Soonpuen  Somapee 

Lieutenant,  Royal  Thai  Navy 

B.S.E.E.,  Royal  Thai  Naval  Academy,  1983 

Submitted  in  partial  fulfillment  of  the 
requirements  for  the  degree  of 

MASTER  OF  SCIENCE  IN  ELECTRICAL  ENGINEERING 

from  the 

NAVAL  POSTGR.'XDUATE  SCHOOL 

September  1988 


ABSTRACT 

This  thesis  describes  the  development  of  computer  algorithms  for  experimental 
measurement  control  and  subsequent  signal  processing  of  transient  signatures  to  syn- 
thesize scattering  impulse  responses  of  scale  model  targets.  The  theories  behind  tran- 
sient scattering  are  considered  in  order  to  construct  the  algorithms.  Up-to-date 
hardware  and  software  technology  are  selectively  implemented  to  optimize  the  resultant 
signal-to-noise  ratio.  The  detailed  explanation  of  the  hardware  and  software  operation 
is  provided.  A  noise  model  for  the  system  is  also  discussed.  Measurement  and  pro- 
gramming validations  are  described,  where  comparisons  are  made  with  numerical  com- 
putations for  transient  scattering  by  selected  canonical  targets.  Results  are  documented 
for  futher  research  in  the  transient  scattering  problem. 
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I.     INTRODUCTION 

A.      BACKGROUND 

In  circuit  theory,  it  is  well  known  that  the  characteristics  of  a  linear  time-invariant 
system  can  be  obtained  by  exciting  the  circuit  with  an  impulse  function.  The  theor\'  of 
transient  scattering  measurements  is  analogous  to  this  method,  and  it  also  provides  a 
unique  parametric  description  of  the  scatterer  via  the  complex  pole  locations  in  the 
Laplace  domain  transfer  function.  There  are  two  approaches  to  obtaining  the  impulse 
response.  One  utilizes  frequency  domain  methods  and  the  other  employs  time  domain 
techniques.  The  former  has  to  sweep  frequencies  over  the  entire  band  but  the  latter  uses 
direct  pulse  excitation.  As  a  result,  the  latter  can  save  time  and  expense.  Therefore, 
direct  transient  scattering  measurements  are  a  viable  alternative  to  the  more  conven- 
tional CW  method  for  obtaining  impulse  responses  of  targets. 

In  1965  Kennaugh  and  iMoffatt  [Ref  1]  investigated  the  properties  of  the  impulse 
response  waveforms  backscattered  from  metallic  targets.  The  approximations  for  tran- 
sient scattering  and  impulse  responses  they  made  formed  the  basis  for  "ramp-response" 
imaging. 

In  1967  DeLorenzo  [Ref  2]  presented  a  system  that  could  be  used  to  obtain  an  ap- 
proximate measurement  of  the  impulse  response  of  scattering  bodies.  His  major  devel- 
opment was  an  indoor  laboratory'  scattering  range  which  facilitated  the  time  domain 
analysis.  He  showed  the  agreement  between  his  calculations  and  the  measurement  of  the 
impulse  response  of  a  rectangular  plate. 

Also  in  1967,  Murchison  and  Falk  presented  a  hybrid  system  for  measurement  of 
impulse  responses  [Ref  3].  This  measurement  system  was  separated  into  two  frequency 
ranges.  For  frequencies  below  4  Ghz,  they  used  an  impulse  generator  while  for  the  fre- 
quencies above  4  Ghz,  they  used  a  swept  frequency  CW  source. 

In  1971  Baum  [Ref  4]  demonstrated  theoretically  that  characteristic  resonances  can 
be  obtained  for  scattering  objects.  His  technique  was  called  the  Singularity  Expansion 
Method  (SEM).  It  was  postulated  that  the  resonances  could  be  obtained  from  the 
transient  scattering  measurement  for  a  complex  object. 

In  1973  Tesche  [Ref  5]  applied  the  SEM  to  to  numerically  evaluate  complex  reso- 
nances of  simple  metallic  scattering  objects. 


In  1975  Blaricum  and  Mittra  [Ref.  6]  proposed  a  technique  for  extracting  the  poles 
and  residues  of  a  system  directly  from  its  transient  response.  Their  technique  provided 
a  number  of  numerical  advantages. 

In  1980  the  original  Transient  Electromagnetic  Scattering  Laboratory  (TESL)  be- 
came operational  at  the  Naval  Postgraduate  School  (NPS).  Hammond  [Ref  7]  devel- 
oped this  laboratory  for  inverse  scattering  research  based  upon  synthesized  ramp 
responses.  This  laboratory  utilized  an  outside  ground-plane  range.  This  range  was 
limited  to  symmetric  scattering  objects  which  could  be  "mirror  imaged". 

In  1983  the  free-field  TESL  began  operation  at  NPS  [Ref  8].  This  laboratory  is  lo- 
cated in  Spanagel  535.  The  TESL  facilitated  research  efforts  in  high  quality  transient 
scattering  measurements  [Ref  9]  used  for  radar  target  identification  and  characterization 
using  complex  pole  natural  resonances  [Ref  10]. 

Prior  to  the  current  undertaking,  McDaniel  [Ref  11]  showed  good  agreement  be- 
tween the  computation  and  the  measurement  of  the  impulse  responses  for  a  sphere  and 
a  thin  wire.  These  measurements  and  signal  processing  were  performed  using  the 
Tektronix  hardware  set  up,  as  originally  employed  in  the  1980  TESL.  In  1987  a  rev- 
olutionary new  Hewlett  Packard  digital  processing  oscilloscope  became  available.  In 
addition,  computational  hardware  had  evolved  dramatically  since  the  early  1980's.  This 
thesis  describes  the  implementation  of  this  new  technology  to  enhance  the  performance 
of  transient  scattering  and  signal  processing  at  the  NPS  TESL. 

B.      OVERVIEW  OF  THIS  EFFORT 

The  main  objective  of  this  thesis  was  to  incorporate  new  measurement  and  compu- 
tational hardware  and  to  design  and  implement  new  software  to  automate  high-quality 
transient  scattering  measurements.  In  this  thesis,  much  work  has  been  done  in  the  de- 
velopment of  the  transient  scattering  work  station  which  has  the  ability  to  measure  and 
synthesize  the  scattering  impulse  response  of  scale  model  targets.  The  major  hardware 
includes  a  microcomputer,  a  digital  programable  oscilloscope  (DPO)  and  a  shielded 
anechoic  chamber.  Automatic  transient  scattering  measurements  have  been  imple- 
mented by  utilizing  an  IEEE-bus  controller. 

With  the  new  HP54120T  DPO,  the  impulse  responses  of  scale  models  are  synthe- 
sized at  a  high  signal-to-noise  ratio  (SNR).  The  acquisition  and  signal  processing  algo- 
rithms run  much  faster  on  the  IBM  PC-AT  vis-a-vis  the  Tektronix  interpreted  BASIC 
computer  previously  used.  The  programs  now  use  QUICK  BASIC  and  FORTRAN 
compilers.     Validations  of  software  were  made  by  way  of  comparing  experimental 


measurements  with  computations  for  simple  canonical  targets.  Various  scale  models 
were  measured  and  doci  mented  in  detail.  Furthermore,  a  noise  calculation  was  provided 
for  further  research. 

This  thesis  is  broken  into  six  chapters.  Chapter  II  describes  system  facilities  and 
hardware.  This  includes  laboratory  description,  anechoic  chamber,  HP54I20T  digital 
processing  oscilloscope  and  the  GURU  II  hardware  interface. 

Chapter  III  presents  the  theory  of  transient  scattering  measurements.  This  is  com- 
posed of  the  free  field  range  system  representation,  mathematical  model  and  transient 
response  solution. 

Chapter  IV  deals  with  acquisition  algorithms  and  signal  processing.  This  chapter 
covers  operational  description,  GURU  II  software  interface  and  the  deconvolution  al- 
gorithm.   In  addition,  some  selected  outputs  are  illustrated. 

Chapter  V  emphasizes  the  experimental  results.  Two  types  of  targets  are  presented. 
One  is  a  canonical  target  and  the  other  is  a  scale  model  target.  The  noise  estimation  of 
the  system  is  also  discussed  in  this  chapter. 

Chapter  VI  summarizes  all  work  and  provides  some  future  considerations. 


II.     SYSTEM  FACILITIES  AND  HARDWARE 

A.  LABORATORY  DESCRIPTION 

Figure  I  on  page  5  shows  the  general  layout  of  the  scattering  range.  It  is  composed 
of 

•  Digital  Programmable  Oscilloscope  (DPO) 

•  Step  pulse  generat'^^r 

•  Low  loss  cables 

•  Broad  band  solid  state  power  amplifier 

•  Broad  band  transmitting  and  receiving  antennas 

•  Anechoic  chamber 

•  Personal  computer  (PC). 

Note  that  the  DPO  includes  the  sampling  and  pulse  generator  unit.  The  DPO  is 
used  as  both  the  master  pulse  source  and  as  the  receiver.  The  antenna  connection  cables 
were  obtained  from  the  DPO  accessory'  kit  for  use  in  low  loss  network  measurements. 
They  are  flexible  for  hard-to-access  areas.  The  characteristics  are  shown  in  Table  21  and 
Table  22  on  page  78  in  Appendix  C. 

Transmitting  and  receiving  antennas  are  of  the  double-ridged  horn  type.  Their  fre- 
quency range  is  1-12.4  GHz,  with  less  than  IdB  of  ripple  in  the  boresight  gain  over  this 
multi-octive  band.  The  Avantek  solid  state  power  amplifier  has  an  advertised  3dB 
bandwidth  of  2  to  6  GHz  but,  by  our  observations,  has  a  usable  bandwidth  of  1  to  7 
GHz.   A  15V  supply  is  needed  for  the  amplifier. 

The  anechoic  chamber  serves  two  purposes.  First,  it  is  used  to  isolate  the  target 
signal  from  the  electromagnetic  noise  and  interference  outside  of  the  chamber.  Secondly, 
the  absorbing  material  on  the  walls,  floor  and  ceiling  of  the  chamber  are  designed  to 
provide  minimal  reflections,  thus  approximating  unbounded  space.  The  computer, 
which  serves  both  acquisition  and  signal  processing  rolls,  is  an  IBM -AT  with  640  K 
RAM,  a  20  megabyte  hard  disk,  and  a  GPIB-PC2A  interface  card.  The  GPIB-PC2A 
was  added  to  the  PC  in  order  to  have  the  abihty  to  talk  or  listen  to  the  HP54120T  DPO. 
The  PC  was  set  to  run  at  9  MHz.  It  has  two  floppy  disk  drives,  and  two  output  ports. 
One  is  a  parallel  port  for  the  printer  and  the  other  is  a  serial  port  for  the  plotter.  The 
color  graphic  board  is  NEC  GBl.   The  display  is  NEC  EGA  color  monitor. 
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Figure  1. 


General  layout  of  the  scattering  range 


B.  ANECHOIC  CHAMBER 

The  anechoic  chamber  is  located  in  room  535  of  Spanagel  Hall.  This  chamber  in- 
corporates a  metallically  shielded  anechoic  enclosure  to  minimize  outside  EM  interfer- 
ence. Chamber  dimensions  are  20  ft.  long  with  a  10  ft.  width  and  10  ft.  height.  Because 
targets  are  suspended  within  the  chamber,  there  is  no  requirement  that  they  have  mirror 
symmetry,  as  is  the  case  for  ground  plane  based  ranges.  In  the  measurements  that  were 
taken  for  this  thesis,  the  target  is  mounted  on  a  thin  styrofoam  column  at  a  distance  of 
8.5  feet  from  the  transmitting  and  receiving  antennas.  Both  receiving  and  transmitting 
antennas  are  placed  close  together  to  closely  approximate  monostatic  measurements. 
The  chamber  interior  is  covered  by  a  special  carbon  impregnated  foam  absorbing  mate- 
rial. The  material  is  cut  into  longitudinal  wedges  along  the  side-walls,  floor  and  ceiling 
to  act  as  conduits  of  EM  energy  towards  the  back  wall.  The  back  wall  is  covered  by  18 
inch  long  pyramids  to  act  as  slow  absorbing  transitions  for  incident  fields.  Further  in- 
formation about  this  free-field  TESL  can  be  found  from  Refs.  9,  11  and  12 

C.  HP54120T  DIGITAL  PROCESSING  OSCILLOSCOPE 

The  Hewlett  Packard  54120T  is  a  fully  programmable  20  GHz,  four  channel, 
digitizing  oscilloscope  with  a  nine  inch  color  display.  It  is  capable  of  automated  meas- 
urements, digital  storage,  and  TDR  measurements.  This  DPO  has  only  been  in  pro- 
duction since  mid- 1987  and  represents  the  state-of-the-art  in  sampling  technology. 
Previous  to  this  thesis  efi^ort,  a  6  year  old  12  GHz  DPO  from  Tektronix  was  used  to 
perform  transient  scattering  measurements  in  the  TESL. 

The  HP-DPO  is  comprised  of  two  major  components.  These  are  the  HP54120A 
Digitizing  Oscilloscope  Mainframe  and  the  HP  54121A  Four  Channel  Test  Set.  It  also 
contains  extensive  self-tests  to  ensure  proper  functioning.  Both  HP  54120A  and  HP 
54121 A  were  calibrated  together  as  a  system  called  the  HP  54120T.  The  DPO  can  be 
remotely  programmed  by  an  HP-IB  (IEEE-488)  controller  interface.  The  HP-IB  com- 
plies with  the  IEEE  488.2  standard.  The  essential  characteristics  of  the  HP  54120T  are 
summarized  in  Table  1  and  Table  2  on  page  7. 

It  should  be  noted  that  the  HP  54 1 21 A  4  Channel  Test  Set  is  sensitive  to 
electrostatic  discharge  (ESD).  Precautions  should  be  taken  to  keep  ESD  from  damaging 
the  equipment.   Details  can  be  found  in  Ref.  13. 

Waveforms  that  are  stored  in  the  DPO  memory  can  be  sent  to  either  an  HP  Thinkjet 
printer  or  an  HP  Plotter,  both  of  which  are  accessed  via  the  HP-IB  bus.  There  are  two 
ways  to  operate  the  instrument.    One  is  manual  operation  via  the  front  panel  and  the 


Table  1.      HP54120T  CHANNEL  CHARACTERISTICS 


Bandwidth 

20, 12.4  GHz 

Number  of  channels 

4 

Full  scale 

8  division 

Maximum  scale 

80  mV/division 

Minimum  scale 

1  mV,'division 

Programable  dc  offset 

±  500  mV 

Input  maximum 

±  2K  dc  +  ac  peak 

Norminal  impedance 

50  ohm 

Connectors 

3.5  mm 

Percent  reflection 

less  than  5%  for  30  ps  risetime 

Table  2.      HP54I20T  TDR  AND  TIMEBASE  CHARACTERISTICS 


Risetime 

less  than  45  psec  at  12.4  GHz  bandwidth  in  average 
mode 

Timebase   full-scale 

10  division 

Minimum  scale  factor 

10  psec'div 

Maximum  scale  factor 

Isec  div 

Time  interval  accuracy 

less  than  10  psec  ±  0.1  %  of  reading 

other  is  to  operate  by  sending  HP-IB  bus  commands  via  the  computer.  A  major  effort 
in  this  thesis  was  developing  programs  to  provide  easy  operator  control  and  signal  ac- 
quisition of  the  DPO  using  the  computer.  These  algorithms  are  explained  in  detail  in 
chapter  IV. 

Ensemble  waveforr^  averaging  has  been  implemented  in  acquiring  scattering  signa- 
tures. The  averaging  reduces  noise  by  \jJW,  where  N  is  the  number  of  averages.  The 
minimum  discernible  signal  level  in  the  digital  processing  oscilloscope  system  is  approx- 
imately 35/1 V.  Typical  noise  amphtudes  for  the  DPO  sampling  unit  are  shown  in  Table 
3  on  page  8. 

Data  display  resolution  is  1024  points  horizontally  x  256  points  vertically  (for  dis- 
play times  of  at  least  200  ps,'div).    The  display  can  be  selected  in  persistence  or  average 


Table  3.      TYPICAL  NOISE 


20  GHz  bandwidth,  Avg  =   1 

Noise  1.2  mV  (RMS) 

20  GHz  bandwidth,  Avg  =  256 

Noise  80/^ F  (RMS) 

12.4  GHz  bandwidth,  Avg  =  1 

Noise  500m K( RMS) 

12.4  GHz  bandwidth,  Avg  =  256 

Noise  35m r  (RMS) 

12.4  GHz  bandwidth,  persistence 

Noise  400m F( RMS) 

mode  and  the  number  of  averages  can  be  specified  as  a  power  of  2,  up  to  2048.  The  data 
output  transfer  rate  and  the  data  record  lengths  are  shown  in  the  following  table. 

Table  4.      DATA  OUTPUT  TRANSFER  RATE  AND  DATA  RECORD  LENGTHS 


10  psec/div  <  time  div  <  20  psecdiv 

100,  or  400  points/record 

20  psec/div  <  time,  div  <  50  psec/div 

100,  400  or  800  points'record 

50  psec  div  <  time,  div  <  200  psec/div 

100,  500  or  1000  points/record 

200  psec  div  <  time/ div  <  1  sec/ div 

128,  256,  500,  512  or  1024  points/record 

Data  output  transfer  rate 

115  kbytes/ s 

For  the  time  domain  reflectometer  (TDR)  system,  the  risetime  is  an  adjustable  value 
based  upon  the  timebase  setting.  The  minimum  is  10  ps  or  .08x  time/div,  whichever  is 
greater.  The  maximum  is  5  x  time  div.  Usually  channel  1  is  used  for  reflection  meas- 
urements, however,  in  this  thesis,  it  was  used  as  the  source  generator. 

The  digitizer  uses  a  12-bit  successive  approximation  A/D  converter.  The  full-scale 
range  of  the  A/D  is  640  mV.  The  least  significant  bit  (LSB)  of  the  A/D  converter  equals 
250  ^V.  This  gives  one  part  in  2560,  or  more  than  11  bits  of  resolution.  Averaging  can 
stretch  out  the  resolution  to  32  ^j.V.   This  improves  the  resolution  to  around  14  bits. 

The  signal  is  sampled  and  digitized  at  a  rate  determined  by  the  trigger  rate,  repe- 
tition rate,  time  base  range,  display  mode,  and  number  of  channels  turned  on.  Further 
information  concerning  the  characteristics  of  the  HP  54120T  can  be  found  from  Ref. 
14. 

D.      GURU  II  HARDWARE 

The  Textronix   supphed  General-purpose  interface  bus  User's   Resource  Utility 
(GURU  II)  is  a  combination  hardware  and  software  package  that  allows  the  use  of  any 


of  the  IBM  personal  computers  (PC),  and  some  PC  compatibles,  to  control  GPIB 
programable  instruments.  In  this  section,  the  focus  is  on  the  hardware.  The  software 
is  explored  in  the  next  chapter. 

The  hardware  consists  of  the  PC2A  GPIB  interface  card  that  plugs  into  the  PC  main 
board,  and  the  GPIB  cable.  The  GPIB  PC2A  interface  card  is  an  8-bit  half-size  board 
which  performs  as  a  talker,  listener  and  controller.  It  also  can  be  implemented  with  the 
full  range  of  talker,  Hstener,  serial  poll,  service  request,  and  remote  programming  func- 
tions for  as  many  as  14  devices  on  the  bus.  It  can  perform  as  a  complete  controller. 
Using  the  PC2A,  data  transfer  between  the  system  memory  and  the  GPIB  are  performed 
at  rates  of  more  than  300  kilobytes  per  second.  The  PC2A  is  compatible  with  all  revision 
levels  of  the  IEEE-488  standard  including  the  HP-IB. 

The  GPIB-PC2A  interface  card  consists  of  these  major  sections: 

•  address  decoding 

•  buffering  and  data  routing 

•  interrupt  arbitration 

•  direct  memory  access  (DMA)  arbitation. 

•  configuration  switches  and  jumpers 

•  GPIB-adapter-TLC  (talker/listener/controller) 

•  time  of  day  clock  with  battery  back-up. 

The  address  decoding  monitors  the  address  lines  to  recognize  when  the  GPIB-PC 
I/O  address  is  present  on  the  computer  I/O  channel  and  enables  read  and  write  access 
to  the  GPIB  adapter  register.  The  buffering  and  data  routing  handles  data  transfer  be- 
tween the  IBM  PC-AT  I/O  channel  and  the  GPIB  adapter  through  a  bidirectional 
internal  data  bus. 

Interrupt  arbitation  recognizes  when  interrupts  have  been  enabled  or  disabled  and 
passes  or  inhibits  them  accordingly.  DMA  arbitation  recognizes  when  DMA  operations 
are  enabled  or  disabled,  and  when  the  last  transfer  has  taken  place.  It  also  routes  the 
DMA  request  and  acknowledges  signals  to  the  selected  DMA  channel.  The  DMA  is 
used  to  transfer  data  directly  between  a  peripheral  and  memory  without  CPU. 

The  IEEE-bus  controller  forms  an  interface  between  the  computer  and  the 
IEEE-488  bus.  Conceptually,  the  IEEE-488  bus  behaves  in  a  way  very  similar  to  the 
system  bus.  Hewlett-Packard  originally  devised  it  [Ref  15]  and  intended  it  to  link  to- 
gether programmable  instruments  in  a  laboratory  or  industrial  environment.  By  con- 
trolling test  equipment  and  measuring  devices  from  the  IEEE-488  bus,  implementing  an 


automatic  testing  station  is  possible.  A  system  under  test  is  connected  to  the  test 
equipment  and  measuring  devices.  The  computer  configures  all  equipment  via  the 
IEEE-488  bus  and  then  reads  the  test  results  fi-om  the  same  bus.  The  GPIB  has  the 
following  limitations: 

•  Half  of  the  devices  attached  to  the  GPIB  must  have  power  on. 

•  The  total  length  of  the  bus  cable  must  not  exceed  20  Meters. 

•  No  more  than  15  devices  may  be  attached  to  one  bus. 

More  detailed  information  can  be  found  in  Ref  16. 
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III.     THEORY  OF  TRANSIENT  SCATTERING  MEASUREMENTS 

A.      FREE-FIELD  RANGE  SYSTEM  REPRESENTATION 

The  physical  model  of  the  system  free  field  scattering  range  was  first  presented  by 
Morgan  in  Ref  12.  The  system  block  diagram  is  shown  in  Figure  2  on  page  12.  Im- 
portant components  and  parameters  of  the  system  are  summarized  in  the  frequency  do- 
main representation  in  the  following  table. 

Table  5.      FREQUENCY  DOMAIN  FUNCTIONS 


Pulse  generator,  driving  the  transmitting  antenna 

XiJ) 

Transmitting  antenna  response 

Hrif) 

Receiving  antenna  response 

H,{f) 

The  multiple  scattering  between  target  and  chamber  wall 

Hsif)  -  -  //c 

Direct  couple  of  receiving  and  transmitting  antenna 

H,{f) 

The  multiple  scattering  between  calibration  sphere  and 
absorber 

Hkif) 

The  multiple  scattering  between  target  and  absorber 

HlciJ) 

The  caUbration  sphere  transfer  function 

m(j) 

The  target  transfer  function 

HKf) 

Thermal  noise 

A^(/) 

The  transient  pulse  source  is  the  step  generator  in  the  TDR  system.  This  generator 
not  only  supplies  the  low-level  transmitted  pulse,  before  amplification,  but  also  acts  to 
coherently  trigger  the  sampling  circuits  of  the  DPO  receiver.  It  thus  is  akin  to  the 
"master  oscillator"  in  a  coherent  pulsed  CW  radar  system.  The  risetime  of  the  step 
generator  is  approximately  25  psec.  This  step  pulse  is  first  amphfied  and  filtered  by  a 
broad-band  GaAs  FET  amplifier. 

The  amplified  pulse,  which  no  longer  resembles  a  step  waveform,  is  then  fed  through 
the  transmitting  antenna  which  radiates  into  the  anechoic  chamber.  The  radiated  field 
emanating  from  the  antenna  is  approximately  an  angular  weighted  spherical  wave.  At 
the  target  position  the  incident  field  is  a  quasi-plane  wave.  Because  of  the  close  prox- 
imity of  the  transmitting  and  receiving  antennas,  some  strong  coupling  of  energy  is  in- 
duced in  the  receiving  antenna. 
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Figure  2.       System  representation  of  the  TESL 
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This  undesirable  signal  can  be  seen  clearly  if  the  time  delay  setting  on  the  DPO  is  ad- 
justed to  begin  many  nanoseconds  before  the  target  response  time-window.  Therefore, 
the  time  setting  delay  must  be  checked  when  the  target  is  present. 

The  initial  cahbration  of  the  beginning  of  the  target  time-window  was  performed  by 
using  a  rectangular  copper  plate,  1  ft  wide  and  1.5  ft  high,  as  the  test  scatterer.  The  back 
scattering  was  observed  with  the  copper  plate  turned  broadside  to  the  centroid  of  the 
horn  positions  and  placed  where  the  target  would  be  mounted.  The  copper  plate  scat- 
tering signature  was  so  large  that  it  could  be  clearly  discerned.  The  results  of  this 
measurement  indicated  a  DPO  time  delay  setting  of  38  nsec  was  needed. 

Multiple  scattering  between  the  target  and  the  chamber  absorber  is  represented  by 
the  two  way  interaction  arrow  connecting  //c(/)  to  //s(/).  The  total  noise  of  the  system 
N{f),  is  due  to  both  thermal  emission  from  the  chamber  and  receiver  noise  in  the  sam- 
pling front  end  of  the  DPO.  These  sources  are  combined  together  and  denoted  as  i\„{/) 
for  the  n-ih  measurement  that  is  made.  To  understand  the  theory  of  transient  scattering 
measurements  adequately,  the  mathematical  model  will  be  presented  in  the  next  section. 

B.      MATHEMATICAL  MODEL 

The  five  key  parameters  in  Table  6  on  page  14  play  a  major  role  in  both  the  physical 
meaning  and  mathematical  sense  of  the  acquisition  and  signal  processing  algorithms. 
These  five  parameters  can  be  related  to  the  system  representation,  described  in  the  pre- 
vious chapter,  by  the  following  three  equations: 

Y^if)  =  Hjif)  //«(/)  m  [HaU)  +  Hcif)  +  Hjij)  +  ifsciD]  +  A',(/)  (3.1) 

>  2(/)  =  /^7(/)  IhiJ)  m  UJa^  +  /^c(/)]  +  A'lC/)  (3.2) 

J  3(/)  =  //7(/)  /^«(/)  M  [Ha^J)  +  //c(/)  +  JlliJ)  +  /4-(/)]  +  A^3(/)-  (3.3) 

1  hese  three  equations  represent  the  three  respective  measurements  that  are  neces- 
sary for  processing  the  scattering  signature  of  a  given  target.  As  shown  below,  these 
measurements  are  for  a  calibration  sphere,  whose  measured  scattering  can  be  compared 
to  an  accurate  calculation;  the  "background"  which  denotes  no  target  present,  and  the 
actual  target  scattering  measurement. 

Equation  3.1  states  that  the  response  to  the  amplified  pulse  generator,  X(0,  due  to 
the  calibration  sphere  and  the  chamber  is  equal  to  the  product  of  the  pulse  excitation, 
transmitting  antenna  transfer  function  H^{f),   receiving  antenna  transfer  function  //«(/) 
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Table  6.      MATHEMATICAL  NOTATION 


'I  he  calibration  measurement  waveform 

y^f) 

The  background  measurement  waveform 

y^if) 

The  target  measurement  waveform 

y^if) 

The  subtracted  cahbration  waveform 

y.(j) 

1  he  subtracted  target  waveform 

K(/) 

times  the  sum  of  the  antenna  coupling  H^{f)  plus  the  calibration  sphere  Hsi/)  ^ind  the 
multipath  multiple  scattering  H§cif)-    Added  to  this  product  is  the  thermal  noise  N„{f). 

Equation  3.2  represents  the  background  measurement,  with  no  target  present.  The 
system  response  is  similar  to  that  previously  considered,  but  with  the  absence  of  the 
transfer  function  for  the  calibration  sphere.  Equation  3.3  expresses  the  target  measure- 
ment in  the  frequency  domain.  The  intepretation  follows  the  same  pattern  as  equation 
3.1  except  the  target  transfer  function,  //Jl/),  appears  instead  of  that  for  the  calibration 
sphere.    The  multipath  multiple  scattering  in  equation  3.3  also  changes  and  is  denoted 

by  //IcOO. 

The  three  equations  as  stated  above  are  three  crucial  measurements  in  the  transient 
scattering  problem.  1  he  backscatter  measurement  with  the  canonical  (calibration)  ob- 
ject }',(/),  measurements  without  target  (background)  Yjif)  and  measurements  with  the 
desired  target  Yiif)  are  obtained  directly  using  the  time  domain  acquisition  program  in 
y\ppendix  B.  This  program  was  written  in  BASIC  as  part  of  this  thesis  elTort.  To  speed 
up  the  execution  of  the  program  by  a  factor  of  about  fifty  times  vis-a-vis  direct  inter- 
pretation, it  has  been  compiled  using  the  Microsoft  Quick  Basic  Compiler. 

After  acquiring  accurate  data,  the  next  step  is  to  process  the  data  and  extract  the 
impulse  response.   This  is  the  topic  of  the  next  section. 

C.      TRANSIENT  RESPONSE  SOLUTION 

The  three  measurements  that  were  just  considered  form  the  basis  for  estimating  the 
"smoothed"  impulse  response  of  the  target.  The  band-limited  nature  of  the  measure- 
ments does  not  allow  synthesis  of  the  true  infinite  bandwidth  impulse  response  which  is 
due  to  an  incident  delta  function  plane  wave.  Instead,  the  "smoothed"  term  means  that 
we  can  estimate  the  transient  target  response  due  to  some  predefined  incident  signal  (a 
double-Gaussian  shaped  waveform  is  used  here)  whose  significant  spectral  content  does 
not  exceed  the  bandwidth  of  the  measurements.    Aside  from  the  currently  used  solid- 
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state  power  amplifier,  the  bandpass  range  of  the  DPO-antenna  system  is  about  1  GHz 
to  12  GHz.  The  spectral  bottleneck  is  due  to  the  GaAS  FET  amplifier,  whose  passband 
is  about  1  GHz  to  7  GHz. 

in  order  to  synthesize  the  smoothed  transfer  fiinction  of  the  target  Hl(f),  the  first 
step  is  to  subtract  the  "clutter"  portion  of  the  responses  in  equations  3.1  and  3.3  which 
is  embodied  in  equation  3.2.   This  subtraction  results  in  equations  3.4  and  3.5  below: 

n(/)  =  Y,{J)  -  \\{J)  (3.4a) 

Y,{J)  =  HjiJ)  Hi,{J)  Xij)  [//loo  +  //fc(/)]  +  A'4(/)  (3.4^) 

A'4(/)  =  A^(/)-A^2(/)  (3.4c) 

J 5(/)  =  ^300  -  >'2(/)  (3.5a) 

Ys^  =  ihU)  liRif)  m  [//]"(/)  +  //icC/)]  +  A'sC/)  (3.5Z;) 

A^5(/)  =  A3(/)  -  A''2(/).  (3.5c) 


Equations  3.4  and  3.5  are  the  very  critical  steps  in  obtaining  the  smoothed  transient 
response  of  the  target.  Ihe  critical  nature  of  the  subtractions  has  to  do  with  coherency 
in  time  of  the  subtracted  clutter.  The  calibration  sphere  measurement  }',(/)  had  the 
background  measurement  }2(/)  subtracted  in  order  to  eliminate  the  clutter.  The  clutter 
signal  is  that  which  is  present  aside  from  the  target  and  is  composed  of  two  main  con- 
stituents. One  is  the  direct  coupling  of  the  receiving  and  transmitting  antennas,  //^(/), 
and  the  other  is  direct  scattering  (as  opposed  to  multipath)  from  the  absorber  in  the 
chamber,  //c(/).  Since  the  measurements  are  made  directly  in  the  time  domain,  the  sub- 
tractions are  done  in  time  domain.     In  a  similar  fashion,  equation  3.5  can  be  obtained. 

The  results,  }4(/)  and  Ys{i),  are  then  transformed  into  frequency  domain  by  discrete 
fast  Fourier  transform  (  FFT  ).  1  his  leads  to  the  next  step,  which  is  to  form  an  optimal 
deconvolution  estimator  for  target  response  given  by 

^oiJ)  ^hin  =  — —- 71 — - —  (3-6) 

\KiJ)\  +c 

where  ]1{J)  is  the  complex  conjugate  of  )\{J),  Yscif)  is  the  computed  calibration  sphere 
and  C  is  the  smoothing  parameter. 
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ysc^  =  Xo{f)Hi{J)  (3.7) 

Xo{t)^FFT--Xoif)  (3.8) 

xoit)  =  Q  e-  *■  ^'-  ''^'  +  C2  ^"  "^^'-  '"^\  (3.9) 

where,  for  a  specified  a,  and  a^  in  this  "Double-Gaussian"  (  DG  )  pulse,  the  C,  and  C^ 
are  chosen  such  that 

^o(0=l  (3.10) 

and 


J. 


+00 

;c^(0^/  =  0.  (3.11) 


Each  of  the  two  Gaussian  waveform  in  equation  3.9  will  be  specified  here  by  a  "10% 
pulse  width",  defined  by  A/^  such  that 


-  at 


^-^*      4     =0.1,  (3.12a) 

which  gives 

a,  =  — ^   ln(10)  =  -f^.  (3.126) 

(A  t^f  (A  t^r 

The  estimator  in  equation  3.6  is  due  to  Riad  [Ref  17]  and  results  from  an  optimal 
least-square  fidelity  criterion  for  deconvolved  target  response.  The  result  of  equation 
3.6,  XJJ)  Hj  ,  represents  the  frequency  domain  scattering  response  of  the  measured  tar- 
get to  an  incident  plane  wave  whose  Fourier  transform  is  XoiJ).  This  incident  waveform 
is  specified  by  numerically  computing  the  scattering  response  of  the  calibration  sphere 
to  it.  This  numerically  generated  response  of  the  calibration  sphere  is  embodied  in 
Yscif)  ,  as  shown  in  equation  3.7.  It  is  obtained  from  the  SPRSCT  program  which  uses 
the  magnitude  and  phase  data  from  the  MIE  series  program.  Both  of  these  programs 
were  written  by  Morgan. 

Both  subtraction  and  optimal  deconvolution  algorithms  are  implemented  in  the 
Deconvolution  Program  in  Appendix  A.  The  procedure  to  operate  this  program  will  be 
explained  in  chapter  IV, 
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IV.     ACQUISITION  AND  SIGNAL  PROCESSING  ALGORITHMS 

A.      GURU  II  INTERFACE 
1.      GPIB  Principal 

This  section  emphasizes  the  GURU  II  software  interface.  In  order  to  under- 
stand the  Acquisition  program  operation,  the  GPIB  principle  has  to  be  understood.  The 
General  Purpose  Interface  Bus  uses  a  twenty-four  conductor  cable.  Sixteen  of  these  are 
used  to  send  the  signal.  Data  is  sent  back  and  fourth  between  the  listener  and  the  talker 
in  a  bi-directional  fashion.  This  means  that  both  transmitted  and  received  data  travel 
on  the  same  line.  Eight  lines  are  used  for  8  bits  of  parallel  data  transfer.  A  byte  is  sent 
by  serial  mode  but  each  bit  is  sent  in  parallel  mode.  Three  lines  act  as  the  coodinator 
between  the  talker  and  the  listener.   These  lines  are  assigned  their  duty  as  follows; 

•  no  data  available  line 

•  data  not  ready  line 

•  data  not  accepted  line. 

The  talker  on  the  bus  controls  the  data  available  on  the  line  and  waits  for  the 
listener's  response  on  the  "not  ready  for  data"  line  or  "data  not  accepted"  line.  The  re- 
maining five  lines  are  the  attention  line,  the  interface  clear  line,  the  remote  enable  line, 
the  service  request  line  and  the  end  or  identify  line.  When  the  attention  line  is  active, 
the  talker  and  the  hstener  will  be  specified.  The  controller  has  full  access  to  five  of  these 
lines  and  has  the  capability  to  directly  address  other  devices  to  command  them  to  talk 
or  listen.  The  interface  clear  line  is  used  to  set  the  turn-on  condition.  When  the  remote 
enable  line  is  asserted,  all  devices  on  the  bus  are  under  computer  control  and  nothing 
external  can  interfere.  The  service  request  line  is  used  by  a  device  to  indicate  a  require- 
ment for  attention  and  to  request  an  interruption  of  the  current  sequence  of  events. 
There  are  two  methods  for  the  controller  to  determine  the  device  and  the  message.  One 
is  the  serial  poll  while  the  other  is  the  parallel  poll.  The  polls  obtain  numerical  codes 
which  show  a  specific  event.  The  serial  poll  is  a  mechanism  that  the  controller  uses  to 
obtain  status  from  individual  devices.  As  each  device  is  polled,  it  returns  a  status  byte 
which  controller  decodes  to  determine  the  device  requesting  service.  The  parallel  poll  is 
a  mechanism  for  accepting   and  decoding  the  general  status  of  as  many  as  eight  devices 
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simultaneously.   The  end  or  identify  line  is  used  by  the  talker  to  indicate  the  end  of  the 
message. 

2.      Basic  Softv^are  Requirements 

There  are  3  files  which  the  GURU  requires  in  order  to  run  the  application  pro- 
gram.    These  files  are 

•  GPIB.COM 

•  IBCONF.EXE 

•  CBIB.OBJ. 

The  primary  GPIB  device  drivers  are  read  in  via  the  CONFIG.SYS  file  at  boot-up  time. 
This  file  is  in  the  root  directory  and  must  include  the  foUov^^ing  DOS  statement: 

•  DEVICE  =  GPIB.COM. 

This  causes  the  GPIB  driver  in  GPIB.COM  to  be  loaded  when  the  system  is  booted. 

The  IBCONF.EXE  is  the  executable  program  to  configure  the  GPIB.COM  file 
containing  the  device  drivers  for  each  instrument.  HPDEV  is  the  name  of  the  device 
driver  for  all  Hewlet-Packard  instruments,  including  the  HP  54120T,  whose  decimal  pri- 
mar\'  address  is  17.  The  IBCONF.EXE  (Interface  Bus  Configure)  program  allows  a 
match  between  the  softwere  and  the  hardware.  The  program  prompts  with  menus  of 
details  about  the  hardware  parameters  associated  with  the  interface  card.  This  device 
for  Hewlett-Packard  is  already  configured  in  the  system.  When  the  IBCONF.EXE  is 
run,  the  device  drivers  are  mapped  for  each  instrument  in  the  display.  By  selecting  the 
HPDEV,  and  editing  it,  the  characteristics  of  the  HPDEV  show  on  the  screen.  This  al- 
lows the  operator  to  change  the  characteristics  to  be  compatible  with  the  hardware. 
After  reconfiguration  of  the  GPIB.COM  file,  the  operator  must  return  to  the  operating 
system  and  then  re-boot  the  PC  for  the  new  parameters  to  be  in  effect. 

CBIB.OBJ  is  the  object  file  which  is  linked  with  any  compiled  BASIC  program 
which  needs  to  interface  with  the  GPIB.COM  device  drivers.  In  our  case,  the  Acquisi- 
tion Program,  written  in  GW-BASIC,  is  first  compiled  using  the  Microsoft  QUICK 
BASIC  compiler.  This  process  creates  an  object  (.OBJ)  file.  The  CBIB.OBJ  object  file 
is  then  linked  with  the  Acquisition  object  file.  The  final  result  is  an  executable  (.EXE) 
file. 

B.      ACQUISITION  ALGORITHM 

The  Acquisition  program  was  written  in  Microsoft  GW-BASIC.    The  hardware  is 
interfaced  with  this  program  by  way  of  the  CBIB  object  file.    This  file  comes  with  the 
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GURU  II  software.  After  compiling  the  source  code,  the  object  file  of  the  Acquisition 
has  to  be  linked  to  the  CBIB  object  file  in  order  to  complete  the  creation  of  the  Acqui- 
sition executable  file.  To  run  the  program,  the  Brun3087  execution  file  is  needed  in  the 
same  DOS  Directory.  This  file  is  supplied  with  Microsoft  QUICK  BASIC.  The  Ac- 
quisition program  uses  a  friendly  and  interactive  menu  format.  The  program  begins  with 
a  displayed  welcome  message,  including  the  time  and  date.  This  screen  also  provides  a 
brief  procedure  on  how  to  measure  the  transient  scattering  response.  The  screen  gives 
operator  the  choice  of  using  the  system  default  or  the  last  response.  The  system  default 
was  designed  for  first  time  operators.  The  last  response  was  designed  for  experienced 
operators.  A  main  menu  is  displayed  on  the  screen  after  the  operator  decides  whether 
to  use  the  system  default  or  the  last  response. 

The  main  menu  contains  many  importance  parameters.  A  first  item  is  the  number 
of  data  points  in  the  acquired  time  series.  This  number  has  to  be  a  power  of  2,  but  not 
more  than  1024.  This  limitation  could  be  modified  in  the  program  by  defining  a  larger 
array  size  up  to  2048  points.  The  second  item  is  the  number  of  sub-averages.  This  pa- 
rameter determines  the  number  of  values  to  average  for  each  time  point  when  the 
HP54120T  is  in  the  average  mode.  This  parameter  has  to  be  an  integer  from  1  to  2048, 
also  in  a  power  of  2.  The  third  item  is  the  number  of  data  blocks.  This  determines  the 
number  of  sets  of  data  in  one  record.  The  fourth  and  the  fifth  items  are  a  user  supplier 
identification  for  the  target  and  the  date  of  measurement.  The  sixth  item  tells  the  oper- 
ator about  the  typical  sampling  time  interval.  The  seventh  parameter  is  the  type  of  the 
waveform  being  acquired  on  the  measurement  being  made.  As  discussed  in  the  previous 
chapter,  there  are  three  scattering  waveforms  which  need  to  be  acquired  for  each  target 
being  considered:  calibration  (sphere),  background  and  actual  target.  The  eighth  item 
tells  the  operator  the  time  window  in  nanoseconds.  The  ninth,  tenth  and  eleventh  items 
are  the  maximum  vertical  setting,  the  data  file  name  output  and  the  auto  time  setting. 
At  the  bottom  of  the  menu,  the  operator  is  asked  if  he  wants  to  change  the  items. 

The  second  menu  appears  on  the  screen  when  the  operator  decides  to  change  some 
item  on  the  main  menu.  This  "overlay"  menu  contains  all  changeable  parameters  and 
allows  a  change  to  these  parameters  interactively.  When  the  changes  are  completed,  the 
operator  can  go  back  to  the  main  menu  by  typing  the  letter  "R".  This  brings  the  oper- 
ator to  the  main  menu  and  displays  the  new  item  on  the  menu.  When  everything  is 
ready,  the  operator  selects  "no  change"  to  begin  to  automate  the  DPO.  The  program 
starts  by  initializing  the  DPO  and  setting  all  specified  parameters  as  directed  by  the  main 
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menu.  The  program  stops  when  the  DPO  has  finished  the  setting  operation.  The  pro- 
gram then  asks  for  the  bandwidth  of  the  measurement.  There  are  two  choices:  20  Ghz 
or  12.4  Ghz.  To  protect  from  errors  in  the  initial  setting,  the  program  stops  for  the  final 
check.  As  soon  as  the  operator  initiates  the  acqisition  cycle,  the  beginning  time  shows 
on  the  screen.  For  a  typical  acquisition  composed  of  2048  averages  and  1  data  block, 
it  takes  about  7  minutes  to  finish  the  operation.  In  the  meantime,  the  computer  waits 
for  the  DPO  to  provide  the  raw  data.  When  the  DPO  finishes  this  operation,  the  ending 
time  is  displayed  and  the  PC  opens  a  temporary  hard-disk  file  for  storing  the  raw  data. 
Before  the  operation  is  ended,  the  header  which  is  to  be  placed  on  the  data  file  is  dis- 
played on  the  screen.  The  operator  can  then  observe  any  error  in  the  data  set,  should 
there  be  any. 

The  header  is  transferred  as  a  separate  file.  The  header  provides  important  charac- 
teristics of  the  data,  such  as  format  type,  number  of  points,  number  of  averages,  x  in- 
crement, X  origin,  x  reference,  y  increment,  y  origin,  y  reference  and  range  setting.  All 
of  these  parameters  allow  the  program  to  recover  meaningful  data  from  the  raw  data. 
The  remainder  of  the  program  constructs  the  standard  data  file  structure  and  the  hard 
copy  output. 

At  the  end  of  the  program,  there  is  a  pause  with  a  query  regarding  the  operator's 
intention  to  continue  with  an  additional  measurement.  To  complete  the  acquisition  cy- 
cle for  one  target,  the  three  measurements  previously  discussed  have  to  be  made  in  order. 
These  three  files  will  then  be  supplied  to  the  deconvolution  program  in  the  next  proc- 
essing step.  For  acquiring  more  than  three  measurements,  the  background  must  be 
measured  for  each  target  so  it  can  be  subtracted  from  the  respective  target  measurement. 
The  background  measurement  occurs  either  immediately  before  or  after  the  target 
measurement.  This  reduces  errors  induced  by  long  term  nonstationary  drifts  in  the  sys- 
tem. The  need  for  such  temporal  coherence  of  the  target-background  measurements  has 
been  verified  experimentally  during  the  course  of  this  work.  Typical  outputs  are  as 
shown  in  Figure  3  on  page  21. 

C.      DECONVOLUTION  ALGORITHM 

The  deconvolution  program  was  written  in  Ryan-iMcFarland  (RM)  FORTRAN  in 
order  to  improve  the  run  time.  The  source  code  is  standard  ANSI  77  FORTRAN  so  it 
may  be  recompiled  under  virtually  any  PC  FORTRAN  compiler  such  as  Microsoft  or 
Lahey.  This  program  was  originally  written  in  Tektronix  BASIC  by  Morgan  [Ref  11] 
and  was  restricted  to  use  on  a  Tektronix  4052  microcomputer. 
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Figure  3.       Deconvolution  input  waveform 

The  use  of  a  standard  FORTRAN-77  compiler  generalizes  the  program  so  it  can  be 
run  on  any  IBM  PC  compatible.  This  algorithm  follows  the  mathematic  deconvolution 
as  described  in  chapter  III.  The  program  has  graphic  capabilities  which  are  supplied  by 
GRAFMATIC  software.  The  procedure  to  compile  the  program  is  described  in  the  RM 
FORTRAN  manual.  This  program  must  be  linked  to  the  GRAFMATIC  library  in  or- 
der to  run  properly. 

The  program  begins  by  displaying  the  standard  input  data  file.  This  data  file  format 
has  been  standardized  for  use  with  all  of  the  software  support  systems.  The  header  of 
each  data  file  contains  the  important  information  such  as  identification,  type  of  meas- 
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urement,  date  of  measurement,  time  of  experiment,  the  maximum  time  scale,  data  block 
and  number  of  points.  Types  of  measurement  are  as  explained  in  the  previous  section. 
The  first  data  file  needed  for  the  deconvolution  program  is  that  for  the  computed  scat- 
tering by  the  calibration  sphere  target.  This  computation  data  is  obtained  from  the  MIC 
series  program  and  SPRCT  program,  written  by  Morgan,  as  described  in  the  previous 
chapter.  The  input  data  required  for  these  programs  are  summarized  in  Table  7  on  page 
23  . 

The  sphere  scattering  computation  is  transformed  into  the  frequency  domain  and 
displayed  on  the  screen  in  order  to  observe  the  bandwidth.  The  bandwidth  of  the  com- 
putation has  to  be  approximately  the  same  as  the  estimated  bandwidth  of  the  measure- 
ment. If  the  bandwidth  is  not  appropriate,  the  operator  should  to  terminate  the 
program  via  a  Ctl-C  and  run  the  MIE  series  followed  by  SPRCT  program  to  provide 
appropriate  data. 

In  the  next  step,  the  program  asks  for  the  calibration  and  background  file  names. 
Both  data  are  read  in  and  plotted  on  the  same  axis  in  order  to  observe  if  time  shifting 
has  occurred,  and  if  so,  to  be  able  to  estimate  how  much.  The  operator  is  then  allowed 
to  shift  one  waveform  either  left  or  right  and  then  re-observe  the  overlayed  waveforms. 
Once  a  satisfactory  alignment  is  observed,  the  calibration  waveform  is  subtracted  from 
the  background  waveform  and  the  difference  waveform  is  transformed  into  the  frequency 
domain. 

The  program  then  stops  and  asks  for  the  target  file  name.  It  then  asks  if  the  oper- 
ator desires  to  use  the  same  background  or  another,  depending  on  which  background 
was  closer  to  the  target  completion  acquiring  time.  In  a  similar  fashion,  the  target  data 
and  the  background  data  are  plotted  on  the  same  axis.  If  there  is  an  observed  drift,  the 
program  allows  the  operator  to  shift  the  target  waveform  to  the  left  or  to  the  right  until 
both  target  and  background  are  lined  up.  Then  the  subtraction  can  be  performed.  The 
difference  waveform  is  then  transformed  into  the  frequency  domain.  Typical  scattering 
data  being  processed  by  deconvolution  program  is  shown  in  Figure  3  on  page  21.  The 
overlay  and  shift  is  shown  in  Figure  4  on  page  23  and  the  subtracted  waveform  is  shown 
in  Figure  5  on  page  24.  The  target  for  this  example  was  a  10  centimeter  long  "thick" 
copper  wire,  and  the  cahbration  target  was  a  standard  0.1025  meter  radius  aluminum 
sphere.  After  the  subtraction,  the  difference  calibration  waveform  and  the  difference 
target  waveform  were  tranformed  into  the  frequency  domain. 
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Figure  4.       Overlay 


Table  7.      CALIBRATION  SPHERE  PARAMETERS 


Sphere  radius  in  meters 

0.1025 

Distance  from  sphere  in  meters 

2.5 

Bistatic  angle  in  degrees 

3.0 

Scattering  plane 

I:  plane 

Time  window  in  nanoseconds 

20.0 

Nyquist  frequency  in  Ghz 

12.775 

No.  frequency  points 

257 
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Figure  5.       Subtraction  waveform 

The  three  frequency  domain  waveforms  were  formed  into  Riad's  optimal  deconvo- 
lution  estimator  given  in  equation  3.6.  The  smoothing  constant  "C"  is  selected  by  the 
operator.  This  constant  is  set  to  optimize  the  denominator  threshold  level  in  the  fre- 
quency domain.  If  it  is  set  too  large,  oversmoothing  will  result  which  reduces  the  spec- 
tral content  in  both  the  low  and  high  frequency  regimes.  If  the  threshold  is  set  to  zero, 
a  "naive"  deconvolution  results.  This  provides  vastly  amplified  noise  in  the  signal  and 
typically  looks  like  "garbage".  A  typical  value  for  C  is  0.1.  Riad's  method  yields  an 
estimator  for  the  scattered  waveform  of  the  target  due  to  an  incident  plane  wave  signal 
having  signal  shape  equal  to  that  which  illuminates  the  computed  calibation  sphere  us- 
ing program  SPRCT.    This  frequency  domain  estimator  is  transformed  to  the  time  do- 
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main  and  the  result  is  displayed.  At  the  end  of  the  program,  the  operator  may  try 
another  smoothing  parameter.  If  the  result  satisfies  the  operator,  it  can  be  stored  in  the 
file  whose  file  name  was  provided  by  the  operator.  The  program  also  allows  the  proc- 
essing of  another  target  directly  without  rereading  the  computed  data  and  the  caHbration 
data. 

D.      EXAMPLE  OUTPUT 

The  output  of  the  deconvolution  program  is  in  both  the  time  domain  and  frequency 
domains,  but  the  following  example  shows  only  the  time  domain  output.  This  output 
can  be  stored  on  a  hard  disk  or  a  floppy  disk. 
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Figure  6.       Thick  wire  measurement 
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The  impulse  response  shown  in  Figure  6  belonged  to  the  signature  of  the  thick  copper 
wire  having  a  90  degree  incident  wave  and  a  90  degree  reflected  wave.  The  number  of 
averages  used  in  the  measurement  was  2048.  The  smoothing  was  C  =  0.1.  There  was 
no  time  drift  observed.  One  data  block  was  used.  Other  parameters  can  be  found  in  the 
following  table. 

Table  8.      90  DEGREE  THICK  WIRE  PARAMETERS 


Wire  length  in  meters 

0.1 

Wire  radius  in  meters 

0.0026 

Incident  angle  in  degrees 

90 

DG  pulse  minimum  10%  spread  in  nanoseconds 

0.28 

DG  pulse  maximum  10%  spread  in  nanoseconds 

0.30 

Time  window  in  nanoseconds 

10.00 

Impact  delay  in  nanoseconds 

2 

Scattering  angle  in  degrees 

90 

Note  that  the  minimum  and  the  maximum  10%  spreads  determine  the  respective  a, 
and  Kj  in  the  DG  waveform  of  equation  3.9. 

First,  the  frequency  spectrum  of  the  smoothed  impulse  response  was  calculated  from 
the  output  and  then  transformed  to  the  time  domain  by  an  FFT  [Ref  18].  The  reason 
for  calculating  the  frequency  response  first  is  that  it  is  more  convenient  to  work  in  the 
frequency  domain  when  using  Raid's  method  to  extract  the  smoothed  impulse  response 
using  equation  3.6 

The  dotted  line  in  Figure  6  on  page  25  is  the  actual  measurement.  Because  the 
measurement  was  contaminated  by  noise,  the  result  did  not  agree  exactly.  The  contin- 
uous line  represents  the  computed  thick  wire  smoothed  impulse  response  due  to  the 
same  excitation.  This  computation  was  performed  directly  in  the  time  domain  using  an 
E-field  integral  equation.  Measurements  for  other  types  of  structures  will  be  illustrated 
in  the  next  chapter. 
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V.     EXPERIMENTAL  RESULTS 

A.      CANONICAL  MODEL  SCATTERING 

In  addition  to  the  10  cm  long  thick  wire  considered  in  Chapter  4,  four  additional 
canonical  targets  were  selected: 

•  10  centimeter  thin  wire  (broadside  and  canted  incidence  angles) 

•  8.1  centimeter  diameter  sphere 

•  12.2  centimeter  diameter  sphere 

•  15.5  centimeter  diameter  sphere 

The  iMie  series  program,  MIE,  was  used  to  compute  the  frequency  domain  scatter- 
ing transfer  function  for  the  various  metallic  spheres.  These  transfer  functions  were  then 
used  as  inputs  to  the  SPRCT  program  to  compute  the  transient  scattering  responses  due 
to  the  specified  Double-Gaussian  (DG)  pulse  incident  plane  wave,  as  specified  by 
equation  3,9.  The  thin  wire  computations  were  performed  using  the  TDIG  program. 
This  program  utilizes  a  numerically  evolved  time-domain  electric  field  integral  equation 
for  the  induced  current  distribution  on  the  wire.  The  scattered  field  is  then  calculated 
by  appropriate  numerical  space-time  integrations  using  this  induced  current.  This  and 
the  sphere  scattering  programs  were  written  by  Morgan.  The  measured  and  computed 
DG  pulse  scattering  responses  are  overlayed  in  each  plot  to  validate  both  the  fidelity  of 
the  experimental  procedure  and  the  accuracy  of  the  deconvolution  algorithm.  The 
characteristics  of  each  canonical  target  scattering  case  are  summarized  in  Tables  9-13 
on  pages  28-33.  A  noteworthy  item  is  that  there  is  a  slight  bistatic  angle  (3  degrees) 
between  incident  and  scattered  aspects.  This  has  been  incorporated  into  the  computa- 
tions. As  can  be  readily  observed,  the  agreement  between  measurements  and  computa- 
tions in  Figures  7-11  on  pages  28-32  is  quite  good  in  all  cases.  The  additive  "rippling" 
of  the  experimental  results  is  due  in  part  to  the  residual  noise  and  clutter  contributions. 
This  latter  pollutant  includes  uncancelled  multipath  interactions  of  the  target  and 
canonical  sphere  scatterers  due  to  imperfect  absorption  by  the  anechoic  chamber  sur- 
faces. Additional  differences  in  the  comparisons  are  due  to  the  approximations  which 
are  innate  to  the  computations. 
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Figure  7.       The  12,2  cm  diameter  sphere  vaUdation. 


Table  9.       12.2  CM  SPHERE  PARAMETERS 


Sphere  radius  in  meters 

0.061 

Distance  from  sphere  in  meters 

2.5 

Bistatic  angle  in  degrees 

3 

Scattering  plane 

E  plane 

Time  window  in  nanoseconds 

20.0 

Narrow  10%  pulse  width  in  nanoseconds 

0.28 

Wide  10%  pulse  width  in  nanoseconds 

0.30 
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Figure  8.       The  10  cm  thin  >vire  vaHdation  at  45  degrees 


Table  10.      45  DEGREE  THIN  WIRE  PARAMETERS 


Wire  length  in  meters 

0.1 

Wire  radius  in  meters 

0.00095 

Incident  angle  in  degrees 

45 

DG  pulse  minimum  10°/'o  spread  in  nanoseconds 

0.28 

DG  pulse  maximum  10%  spread  in  nanoseconds 

0.30 

Scattering  angle  in  degrees 

135 

No.  wire  segments 

20 
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Figure  9.       The  10  cm  thin  wire  vaUdation  at  90  degrees 


Table  11.      90  DEGREE  THIN  WIRE  PARAMETERS 


Wire  length  in  meters 

0.1 

Wire  radius  in  meters 

0.00095 

Incident  angle  in  degrees 

90 

DG  pulse  minimum  10%  spread  in  nanoseconds 

0.28 

DG  pulse  maximum  10%  spread  in  nanoseconds 

0.30 

No.  wire  segments 

20 

Scattering  angle  in  degrees 

90 
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Figure   10.       The  8.1  cm  diameter  sphere  validation 


Table  12.      8.1  CM  SPHERE  PARAMETERS 


Sphere  radius  in  meters 

0.0405 

Distance  from  sphere  in  meters 

2.5 

Bistatic  angle  in  degrees 

3.0 

Time  window  in  nanoseconds 

20.0 

Nyquist  frequency  in  Ghz 

12.775 

Narrow  10%  pulse  width  in  nanoseconds 

0.28 

Wide  10%  pulse  width  in  nanoseconds 

0.30 
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Figure   II.       The  15.5  cm  diameter  sphere  validation 


B.      SCALE  MODEL  SCATTERING 

The  scattering  responses  of  the  wire  and  sphere  targets  due  to  the  incident  DG  in- 
cident pulse  can  be  numerically  evaluated  using  straightforward  algorithms.  This  is  not 
the  case  for  more  complex  scattering  objects,  such  as  a  highly  realistic  model  of  an  air- 
craft. Primitive  computational  algorithms  do  currently  exist  for  such  complex  structures 
but  the  poor  accuracy  of  their  results  would  be  insufficient  to  properly  validate  the 
measurements  being  made  here.  On  the  other  hand,  the  very  good  comparisons  shown 
in  the  previous  section  for  both  low-damped  (wires)  and  high-damped  (spheres)  strongly 
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Table  13.       15.5  CM  SPHERE  PARAMETERS. 


Sphere  radius  in  meters 

0.0775 

Distance  from  sphere  in  meters 

2.5 

Bistatic  angle  in  degrees 

3.0 

Scattering  plane 

E  plane 

Time  window  in  nanoseconds 

20.0 

Nyquist  frequency  in  Ghz 

12.775 

No.  frequency  points 

257 

support  the  conjecture  that  high  quality  measurements  can  be  obtained  for  all  metallic 
targets  having  comparable  physical  dimensions. 

After  the  validations  were  made  using  wires  and  spheres,  extensive  measurements 
were  made  of  the  transient  backscaitering  from  four  different  silver  coated  1,72  scale 
plastic  aircraft  models.  The  full-size  dimensions  of  each  aircraft,  whose  identities  will 
not  be  given  here,  are  shown  in  Tables  14-17  on  pages  34-37. 

Measurements  were  made  for  various  incident  aspects  in  the  wing-fuselage  plane  on 
each  model,  using  horizontal  incident  polarization:  0  degrees  (nose-on);  30  degrees;  90 
degrees  (broadside);  and  180  degrees  (tail-on).  In  addition,  a  broadside  look-down  as- 
pect was  considered,  with  backscattering  measured  for  incident  linear  polarizations  both 
parallel  to  ("vertical  wing")  and  perpendicular  to  ("horizontal  wing")  the  fuselage. 
Samples  of  these  numerous  measurements  are  displayed  in  Figures  12-15  on  pages  34-37. 

The  complete  library  of  waveforms  is  being  used  to  support  research  in  demon- 
strating the  accuracy  and  aspect  independence  of  natural  resonance  radar  target  iden- 
tification. 

Four  scale  model  aircraft  were  selected.   They  are  as  follows: 

•  Scale  model  1 

•  Scale  model  2 

•  Scale  model  3 

•  Scale  model  4. 

C.      SYSTEM  NOISE  ESTIMATION 

Noise  in  the  data  acquired  from  the  acquisition  program  was  investigated  in  this 
effort  in  order  to  develop  a  mechanism  to  estimate  the  signal-to-noise  ratio  (SNR).  The 
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Figure   12.       Vertical  and   horizontal  wing  scale  aircraft  model  1 


Table  14.      DIMENSIONS  OF  SCALE  AIRCRAFT  MODEL  1 


Scale 

1/72 

Length 

36  feet  4  inches 

Wingspan 

31  feet 
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Figure  13.       Scale  aircraft  model  1  at  000,  030,  090,  180  degree  aspect  angle 


Table  15.      DIMENSIONS  OF  SCALE  AIRCRAFT  MODEL  2 


Scale 

1/72 

Overall  length 

53  feet 

Extended  wingspan 

28  feet  7  inches 
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Figure   14.       Scale  aircraft  model  comparisons  at  000  degree  aspect  angle 


Table  16.      DIMENSIONS  OF  SCALE  AIRCRAFT  MODEL  3 


Scale 

1/72 

Length 

56  feet 

Wingspan 

40  feet  8  inches 

Height 

15  feet  4  inches 
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Figure  15.       Scale  aircraft  model  comparisons  090  degree  aspect  angle 


Table  17.      DIMENSIONS  OF  SCALE  AIRCRAFT  MODEL  4 


Scale 

1/72 

Length 

63.8  feet 

Wingspan 

42.8  feet 

Height 

18.6  feet 
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SNR  is  a  major  determining  factor  for  the  measurement  fidelity.  Much  of  this  SNR 
assessment  involved  the  computing  some  statistics  of  the  waveform.  All  of  these  statis- 
tical calculations  were  performed  using  GIIAFSTAT  with  APL. 

In  the  acquisition  step,  the  calibration,  the  background,  the  target  and  a  second 
background  were  measured  consecutively.  The  two  backgrounds  were  measured  for  the 
purpose  of  noise  estimation.  In  the  following  analysis,  a  2048  ensemble-average  mode 
was  used  to  acquire  the  data.  The  acquisition  of  each  waveform  took  about  9  minutes 
to  complete. 

In  a  separate  investigation,  the  stationarity  of  the  system  noise  was  investigated  by 
acquiring  seven  consecutive  background  measurements.  These  backgrounds  were  then 
subtracted  from  each  other  to  observe  the  noise  occuring  during  the  consecutive  meas- 
urement. For  example,  the  second  background  and  third  background  were  subtracted, 
as  is  shown  in  Figure  16  on  page  39  along  with  various  statistical  parameters  in  Table 
18  on  page  39  mentioned,  the  statistics  of  each  subtracted  waveform  was  calculated  by 
the  GRAFSTAT  software  packgage.  The  variance  of  the  noise  can  be  calculated  from 
the  standard  deviation  and  expected  power  can  be  approximated  by  this  variance  when 
the  noise  has  a  small  mean.  This  approximation  was  reasonable  because  the  squares  of 
the  mean  for  the  actual  measurement  data  are  very  small,  (  see  Figure  16  on  page  39  ). 
The  computed  noise  power  is  plotted  versus  time  in  Figure  17  on  page  40.  As  can  be 
seen,  the  measured  noise  power  varied  only  slightly  with  time.  A  linear  least-squares  fit 
yields  the  equation 

Y=  -93.042  +  0.035727A'  ,  (5.1) 

where  Y  is  the  expected  noise  power  and  X  is  the  time  in  minutes.  This  model  was 
based  on  only  seven  background  measurements  with  only  2048  averages  per  measure- 
ment. For  practical  purposes  the  system  noise  can  be  assumed  to  be  fixed  in  time,  de- 
pending only  upon  the  number  of  ensemble  averages  used  to  complete  the  waveform 
measurement. 

To  calculate  the  SNR,  the  thick  wire  measurement  was  used.  The  signal  plus  noise 
waveform  is  displayed  in  Figure  18  on  page  41  while  the  computed  statistical  parameters 
are  given  in  Table  20  on  page  41. 

This  waveform  was  time  windowed  to  eliniinate  the  noise  outside  of  the  window 
before  the  statistical  parameters  were  calculated.  The  signal  plus  noise  power  was  cal- 
culated by  way  of  the  variance  using  the  following  equation. 
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Figure   16.       Noise  waveform  statistics 


Table   18.      STATISTICS  OF  NOISE  WAVEFORM 


Mean 

-3.0211E-10 

Std.  deviation 

0.000023345 

Median 

-7.019E-8 

Maximum 

0.000066731 
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Figure  17.       Noise  power 


Table  19.      STATISTICS  OF  NOISE  MODEL 


Mean 

-91.917 

Std.  deviation 

0.99986 

Median 

-92.343 

Maximum 

-90.669 
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Figure   18.       Signal  plus  noise  waveform  statistics 


Table  20.         STATISTICS     OF     THICK     WIRE     SIGNAL     PLUS     NOISE 
WAVEFORM 


Mean 

-6.2991E-6 

Std.  deviation 

0.0003362 

Median 

0 

Maximum 

0.0014128 
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E{X^)  =  Variance  +  E\X)  (5.2) 

where  E{X)  is  mean  of  the  noise  and  E{X^)  is  the  mean  square  value  of  the  noise  (noise 
power).  Note  that  the  second  term  on  the  right  side  of  equation  5.2  was  close  to  zero, 
so  it  can  be  neglected.  As  a  result,  the  expected  power  is  equal  to  the  variance.  The 
signal  plus  noise  power  was  equal  to  0.0003382^  or  -69.41  dB.  The  noise  expected  power 
was  calculated  by  equation  5.1  and  is  given  by  -92.64  dB  for  the  average  time  of  10 
minutes.  Since  the  noise  power  was  so  small  compared  to  the  signal  plus  noise  power, 
the  subtraction  between  the  signal  plus  noise  and  the  noise  was  not  necessary.  The 
noise  power  when  time  windowed  was  reduced  by  a  factor  of  200/512,  or  4.08  dB.  The 
result  for  the  SNR  is  27.31  dB  in  the  case  of  the  thick  wire  measurement.  Note  that  this 
SNR  is  about  10  dB  better  than  that  provided  for  the  same  target  using  the  Tektronix 
equipment,  as  carried  out  by  McDaniel  [Ref  11]. 
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VI.      CONCLUSIONS 

A.      SUMMARY 

This  thesis  describes  the  theory  and  implementation  of  a  newly  updated  transient 
scattering  measurement  facility.  An  overview  of  the  historical  precedents  and  moti- 
vations for  this  development  was  given  in  Chapter  I.  This  was  followed  in  Chapter  II 
by  a  detailed  description  of  the  new  equipment  being  used  as  well  as  the  physical  struc- 
ture of  the  anechoic  chamber. 

Building  upon  this  overview  of  the  hardware  aspects,  the  conceptual  and  math- 
ematical modeling  of  the  scattering  range  is  undertaken  in  Chapter  III.  Using  the 
physical  model  of  the  radiating  and  receiving  antennas  in  the  anechoic  chamber,  the 
various  electromagnetic  interactions  between  the  scattering  target  and  the  chamber  are 
represented  by  a  comprehensive  linear  system  model.  The  topology  of  the  various  fre- 
quency domain  transfer  functions  in  this  model  is  then  used  to  justify  the  measurement 
and  pre-processing  steps  required  to  estimate  the  smoothed  scattering  impulse  response, 
due  to  a  specified  double-Gaussian  incident  plane  wave.  There  are  three  measurements 
required  for  each  target  impulse  response  estimation.  Aside  from  the  actual  target 
measurement,  a  background  measurement  is  needed  to  subtract  out  the  directlv  received 
signal  from  the  transmitting  antenna,  in  addition  to  that  scattered  from  the  various 
chamber  surfaces.  A  third  measurement,  for  a  metallic  sphere,  is  needed  to  allow  for 
compensation  of  the  transfer  function  characteristics  of  the  transmitting  and  receiving 
antennas  in  the  chamber.  By  effecting  a  deconvolution  comparison  with  a  highly  accu- 
rate numerical  calculation  of  the  transient  scattering  response  of  the  sphere  target,  this 
fmal  measurement  also  permits  the  estimation  of  the  target's  response  to  the  same 
excitation  as  is  used  in  the  sphere  computation. 

Chapter  IV  considers  the  details  of  the  actual  software  that  performs  the  signal  ac- 
quisition, signal  averaging  and  deconvolution  post-processing  of  the  measured 
waveforms.  The  acquisition  and  signal  averaging  software  is  a  blend  of  a  commercial 
IEEE-488  bus  controller  (GURU-II)  and  a  custom  driver,  designed  as  part  of  this  thesis 
effort.  The  deconvolution  algorithm  provides,  first  of  all,  graphically  displayed  iterative 
background  subtractions  from  both  the  target  and  sphere  measurements  to  allow  com- 
pensation for  time  drifts.  Secondly,  the  final  deconvolution  is  performed  via  a  variant 
on  Riad's  method,  as  discussed  in  Chapter  III. 
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The  validations  shown  in  Chapter  V  speak  for  themselves  as  to  the  expected  meas- 
urement quality  obtained  using  this  new  transient  scattering  range.  Both  low-Q  and 
high-Q  type  targets  are  considered  by  the  respective  comparisons  of  sphere  and  thin-wire 
scattering  measurements  with  computed  smoothed  impulse  responses.  Further  illus- 
trations of  the  type  of  measurements  to  be  performed  with  this  range  are  shown  in  the 
scale  model  aircraft  backscattering  signatures.  A  simple  software  modification  will  allow 
the  display  of  very  broadband  frequency  domain  radar  cross  sections  (RCS)  for  scale 
model  targets. 

To  quantitatively  investigate  the  fidelity  of  the  measurement  system,  a  noise  esti- 
mation of  the  new  implementation  was  investigated  and  was  modelled.  The  particular 
signal-to-noise  ratio  (S\R)  of  a  thick  wire  scattering  measurement  was  calculated.  A 
27  dB  SNR  of  this  case  demonstrated  the  relatively  high  fidelity  of  the  measurement 
system  for  the  case  small  RCS  targets,  including  silver  coated  1/72  scale  model  aircraft. 

B.      FUTURE  CONSIDERATIONS 

Even  though  the  demonstrated  system  performance  was  excellent,  there  is  no  such 
thing  as  too  high  an  SNR  in  scattering  measurements.  In  this  regard,  the  power- 
bandwidth  product  has  the  potential  to  be  enhanced  by  future  efforts. 

The  1  to  7  GHz  passband  of  the  scattering  system  was  limited  by  the  transmitting 
GaAS  power  amplifier.  On  the  other  hand,  the  bandwidth  of  the  HP  DPO  system  can 
go  as  high  as  20  Ghz,  while  the  current  antennas  have  a  1  to  12  GHz  passband.  Future 
use  of  broader  bandwidth  antennas  and  power  amplifier,  having  a  rated  RMS  output 
power  of  2  watts  or  more,  would  increase  the  overall  system  performance.  More  accu- 
rate smoothed  impulse  responses  could  be  achieved,  as  well  as  wider  bandwidth  RCS 
measurements. 

Because  of  the  limited  scope  of  this  thesis,  the  noise  power  estimation  was  only  an 
approximation.  A  more  comprehensive  noise  analysis  should  be  carried  out  in  order  to 
quantify  and  minimize  the  noise  power  sources. 

The  final  consideration  is  the  future  investigation  of  alternate  mechanisms  for  pro- 
viding broadband  scattering  measurements.  One  such  mature  technique  is  through  a 
stepped  frequency  continuous  wave  system,  employing  a  network  analyzer  front-end  as 
the  coherent  receiver,  perhaps  preceded  by  a  low-noise  amplifier  to  reduce  the  overall 
system  noise  figure.  A  somewhat  radical  technique,  one  that  is  currently  being  investi- 
gated at  NFS,  is  to  use  a  very  broadband  amplified  noise  source  as  the  transmitting 
source,  while  employing  the  HP  DPO  as  a  two-channel  correlation  receiver.     This 
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method  is  akin  to  a  dynamic  matched  filter,  ofTering  enhanced  bandwidth,  at  relatively 
low  cost  vis-a-vis  the  stepped  fi^equency  setup.  Additionally,  the  potential  exists  for  both 
tactical  and  strategic  employment  of  such  random  source  systems,  either  as  active 
"pseudo-jammers"  or  as  totally  passive  radars  relying  upon  either  natural  or  man-made 
random  signal  sources  for  the  illumination  of  targets. 
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APPENDIX  A.      THE  IMPULSE  RESPONSE  DECONVOLUTION 

The  following  computer  program  calculates  the  Impulse  Response  from  the  meas- 
urement data  as  described  in  chapter  4.   The  program  is  written  in  RM/ FORTRAN. 

The  subroutine  PL08  and  Pll  require  the  GRAFMATIC  routine  library.  The  pro- 
cedure to  compile,  link  and  run  this  source  code  refers  to  the  RM/FORTRAN  User's 
Guide  Version  2.11  (DOS). 

C  THIS  IS  THE  IMPULSE  RESPONSE  DECONVOLUTION  * 

C  WRITTEN  BY  LT. JG  SOONPUEN  SOMAPEE  Feb  88  * 

C  ADAPTATION  OF  ORIGINAL  TEKTRONIX  PROGRAM   * 

C  BY  PROF.  M.A.MORGAN  Circa  1983-85         * 


c 

WHERE: 

* 

c 

DATE 

IS 

THE  DATE  OF  MEASUREMENT 

* 

c 

DATAB 

IS 

THE  NUMBER  OF  DATA  BLOCKS 

* 

c 

Fl 

IS 

THE  DUMMY  STRING 

* 

c 

ID 

IS 

THE  IDENTIFICATION  STRING 

* 

c 

K 

IS 

THE  NUMBER  OF  POINTS 

* 

c 

L9 

IS 

THE  RELATIVE  SMOOTHING 

* 

c 

MAG 

IS 

THE  MAGNITUDE  OF  Xo*Ht(f) 

* 

c 

NAME4 

IS 

THE  NAME  OF  THE  TARGET  DATA 

* 

c 

N 

IS 

THE  NUMBER  OF  POINTS 

* 

c 

NUMPT 

IS 

THE  SAME  AS  K 

* 

c 

NEXPT 

IS 

INDICATEES  POWER-OF-2  EXPONENT 

* 

c 

WILL  BE  GENERATE  BY  PREFFT 

* 

c 

SET  TO  -1  INDICATES  ERROR 

* 

c 

R 

IS 

THE  REAL  PART  OF  THE  INVERSE 

* 

c 

FOURIER  TRANSFORM  OF  Xo*Ht( f ) 

* 

c 

SUBAV 

IS 

THE  NUMBER  OF  SUB -AVERAGES 

* 

c 

T 

IS 

THE  TIME  INTERVAL 

* 

c 

TIME 

IS 

THE  TIME  INTEVAL 

* 

c 

TITLE 

IS 

THE  TITLE  NAME  FOR  PLOTTING 

* 

c 

TUT 

IS 

THE  FREQUENCY  DOMAIN  OF  Xo*H(f)* 

c 

TYPE 

IS 

THE  WAVEFORM  TYPE 

* 

c 

W 

IS 

THE  COMPLEX  EXPONENTIAL  ARRAY 

* 

c 

GENERATED  BY  PREFFT  SUBROUTINE 

* 

c 

Xo*Ht(f) 

IS 

THE  GAUSSIAN  WEIGHTED  COSINE 

* 

c 

XS 

IS 

THE  STRING  FOR  PLOTTING  ROUTINE* 

c 

FOR  THE  HEADER 

* 

c 

DRIVER 

* 

c 

Y4 

IS 

T;^E  CALIBRATION  SPHERE 

* 

c 
c 
c 

SUBTRACTED  WAVEFORM  (FREQUENCY) 

* 
it 

Y5 

IS 

THE  TARGET  SUBTRACTED  WAVEFORM 

■k 

c 
c 
c 

FREQUENCY  DOMAIN 

* 

•A- 

YSC 

IS 

THE  COMPUTED  SCATTERED  ELECTRIC* 
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c 

FIELD 

c 

Yl 

IS  THE 

c 

Y2 

IS  THE 

c 

Y3 

IS  THE 

c 

YS 

IS  THE 

c 

YN 

IS  THE 

c 

XMIl 

IS  THE 

c 

XMAl 

IS  THE 

FOR  THE  CANONICAL  TARGET  * 
CALIBRATION  MEASUREMENT  * 
BACKGROUND  MEASUREMENT  * 
TARGET  MEASUREMENT  * 
SAME  AS  XS  * 

YES  OR  NO  STRING  * 

MINIMUM  VALUE  * 

MAXIMUM  VALUE  * 

CHARACTER'>64  TITLE ,  XS ,  ID ,  YS ,  DATE ,  TYPE ,  NAME 

CHARACTER*64  TYPEl ,TYPE2,TYPE3,F1 

CHARACTER*16  NAME 1 , NAME 2, NAME 3 

CHARACTER*64  ID1,ID2,ID3 

CHARACTER* 1  YNS , YNS 1 

COMPLEX  YSC( 1025) ,Y4( 1025) ,Y5( 1025) ,W( 1025) 

REAL  T,T2  ,TIME 

REAL  Yl( 1024), Y2( 1024), Y3( 1024) 

INTEGER  N,SUBAV,DATAB,NUMPT 

YNS='N' 

C     READ  DATA  FOR  COMPUTED  SPHERE  THEN 

C     CALLING  SUBROUTINE  FFF  TO  DRAW  THE  TIME 

C     SERIES  AND  SPECTRAL  PLOT 


C**Vc?V*'jViV-V*yf?VVfVc?VVr*5V-V*iViV>ViV**VrVrVf>V**iV**?V***Vf*5ViV**'jV**** 

WRITE('^*) 

'***y«ViV*V«ViV******Vf*Vc**********Vf>V*******iV**TV*Vf)V***TV-jV**' 

VRITE(*,*) 

Sv 

*• 

WRITE(*,*) 

Sv 

DECONVOLUTION  PROGRAM  BY  DR.  M.  A.  MORGAN 

Vf' 

WRITEC'V,*) 

Sv 

ADAPTED  TO  IBM  BY  LT.  JG.  SOONPUEN  SOMAPEE 

,v' 

WRITE(*,*) 

WRITE(*,*) 

Sv 

FORMAT  OF  DATA  FILE: 

*' 

WRITE(*,*) 

WRITE(*,''0 

•* 

1. 

READ  ID              (FORMAT  CHARACTER*64) 

*• 

WRITE (*,*) 

•* 

2. 

READ  TYPE            (FORMAT  CHARACTER''^64) 

*' 

WRITE  ('•%*) 

'yr 

3. 

READ  DATE            (FORMAT  CHARACTER''^64) 

*• 

WRITE(*,*) 

'* 

4. 

READ  Fl              (FORMAT  CHARACTER'V64) 

,v' 

WRITE  C*,'0 

'* 

5. 

READ  SUB -AVERAGE      (FORMAT  INTEGER) 

Vf' 

WRITE (*,*) 

•* 

6. 

READ  TIME  WINDOW      (FORMAT  REAL(F5.2)) 

*' 

WRITEC*,*) 

'Vf 

7. 

READ  DATA  BLOCK       (FORMAT  INTEGER) 

Vf' 

WRITE (*,*) 

'* 

8. 

READ  NUMBER  OF  POINTS  (FORMAT  INTEGER) 

*' 

WRITE (*,*) 

'  yfV« 

f^fsW? 

yfVfVryfyjyryfy-yryfVfVryfyfy-yrVfyfV-yrVfVfyrV-VfTiVyfV-Vry-yfyfyrVfVryryrVryfycyfyfyfVr' 

V*' 

WRITE(''S*) 

'FILE  NAME  FOR  COMPUTED  CALIBRATION  SPHERE:   ' 

READ(*,173) 

NAME 

0 

c 

C 

c 
c 
c 

C 

YSC  =  COMPUTED 

SPHERE  COMPLEX  FREQUENCY  DOMAIN  DATA 

NUMPT  =  NUMBER 

OF 

POINTS  IN  THE  OUTPUT  OF  THIS  ROUTINE 

TIME  =  OUTPUT  TIME 

INTERVAL 

CALL  FFF 1( NAME, YSC, W,NEXPT, TIME) 
9110  CONTINUE 
C     READ  THE  DATA  FILES  FOR  CALIBRATE  AND  BACKGROUND 
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C  THEN  OVERLAY  THESE  ONTO  THE  SAME  AXIS. 

C  PLOT77  ALLOWS  THE  OPERATOR  TO  SHIFT  ONE  WAVEFORM 

C  PRIOR  TO  SUBTRACTING  THE  BACKGROUND. 

C 

C  THE  RESULT  FROM  SUBTRACTION  IS  TRANSFORMED  TO 

C  FREQUENCY  DOMAIN  AND  SAVED  TO  THE  COMPLEX  ARRAY  Y4 

C 

C  NAMEl  =  NAME  OF  THE  FILE  TO  SENT  TO  PLOT77  ROUTINE 

C  (CALIBRATE  DATA) 

C  NAME2  =  NAME  OF  THE  BACKGROUND  DATA 

C  IDl   =  IDENTIFICATION  CHR$  FOR  CALIBRATE 

C  Y4    =  COMPLEX  ARRAY  OUTPUT  (CALIBRATE  ) 

WR I TE  (  *  *  )  '  ***iV Vr?V**?V?V**VMV****^V*TV**iVi^**iV*iViV*****TiV****  ' 

WRITE(*i*)  'FILE  NAME  FOR  CALIBRATION  SPHERE:   ' 

READ(*,173)  NAMEl 

CALL  INPP(NAME1, IDl, TYPEl, DATE, F1,SUBAV,T2,DATAB,NUMPT,Y1) 

WR  I  TE  (  *  *  )   '  iV'5V?V»V*'^^V?V**Vr^ViV**iV**Vf  ?V*yrTV^ViV-.V'A'*:iVVf  ** •>>*** ^VSViV *•>> ' 

WRITE(*!*)  'FILE  NAME  FOR  SPHERE  BACKGROUND  WAVEFORM:  ' 
READ(*,173)  NAME2 

CALL  INPP(NAME2,ID2,TYPE2,DATE,F1,SUBAV,T2,DATAB,NUMPT,Y2) 
CALL  PL077(Y1,Y2,T2,NUMPT,ID1) 

9112  CONTINUE 

TOiTE(*%V)  'FILE  NAME  FOR  TARGET  WAVEFORM:     ' 
READ(*,173)  NAME3 

CALL  INPP(NAME3,ID3,TYPE3,DATE,F1,SUBAV,T2,DATAB,NUMPT,Y3) 

WRITE(*,''0 

WRITE(*,''0  'USE  PREVIOUS  BACKGROUND  FOR  THE  TARGET  ?  (Y/N):    ' 

READ(--,173)  YNSl 

IF(YNS1.EQ.  'Y'.OR.  YNSl.EQ.  'y')  GO  TO  9222 

WRITE (^'S^O  'FILE  NAME  FOR  TARGET  BACKGROUND  WAVEFORM:  ' 
READ('S173)  NAME2 

CALL  INPP(NAME2,ID2,TYPE2,DATE,F1,SUBAV,T2,DATAB,NUMPT,Y2) 

9222  CONTINUE 
C 

C  REPEATING  SHIFT  AND  SUBTRACT  OPERATIONS  ON  TARGET  AND  BACKGROUND 

C  NAME1=  TARGET  DATA  FILE 

C  NAME 2=  BACKGROUND  DATA  FILE 

C  OUTPUT 

C  IDl  =  IDENTIFICATION  OF  THE  TARGET 

C  Y5  =  COMPLEX  ARRAY  OF  FREQUENCY  DOMAIN  TARGET  DATA 

CALL  PL077(Y3,Y2,T2,NUMPT,ID3) 
DO  987  I=1,NUMPT 

Y4(I)=CMPLX(Y1(I),0.0) 
Y5(I)=CMPLX(Y3(I),0.  0) 
987        CONTINUE 
MODE=0 
CALL  FFT( NUMPT , MODE , TIME , NEXPT , W , Y4 ) 
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CALL  FFT( NUMPT , MODE , TIME , NEXPT , W , Y5 ) 

WRITE  f'"''  *^   '  •'f's'f  "■''''i'f*"''i^<'<''!'f'5'f**"**'''f'''f''f*''f**"****'''f**''^f****'^f**** 

CALL  DEC0N3(NUMPT,YSC,Y4,Y5, IDS, TYPES, DATE, F1,SUBAV,T2,DATAB) 

WRITEC*!*)  'WOULD  YOU  LIKE  TO  PROCESS  ANOTHER  TARGET  ?  (Y/N): 

READ(''S173)  YNS 

IFCYNS.EQ.  'Y'  .OR.  YNS.EQ.  'y'  )  GOTO  9112 

WR I TE  (  *  *  )   '  ***if*icicifkiclficific-kieicific-ic-kicicici--k*icirirk-k*  ' 

WRITE(*,*)   '^  eeeee    eeeeee     e 

'*   eee    eee    eeee  *' 

eeeeeeeeeeeeeeeeeeeeeeeeeeeee  *' 

eeeeeeeeeeeeeeeeeeeeeeee  ** 

eeeeeeeeeeeeeeeeeeeeee  *' 

^J^  «M  ^^  mS^  mj^  W«  k.'d  «J^  «.'«  ^«  «Jd  kJ^  hJ^  W^  w^  ^^  k'tf  W^  ^^  ^^  ^d  J^  tJ^  «J^  «J^  W^  h7^  mJ^  mJ^  ^^  *J^  W^  w^  ' 


WRITE(*,^0 
WRITE(*,*) 
WRITE(*,''0 
WRITE (*,*) 

WRITE(*,^0 
WRITE  C*,'^) 
173   FORMAT(A) 
STOP 
END 

SUBROUTINE  INPP(  NAME,ID,TYPE,DATE,F1 ,SUBAV,T2,DATAB ,NUMPT,Y) 

I     TTTTTrTr  **  «k  ««  ««  *»  *K  **1\  *\  #f7*7*  **  #C^«C^«v  *CT*  **  #«  #v  *r<*  **  *v  «?<r^v  7*  *»  7v  7v  7r7T  #v  #^7*  «»  7*  *v  *r7»  *»  «\7i«v  *f  7V  *»  7»  7r«f  *\  #\  7v  7v  #*  *S**  «v  ««  *t7s  7T7* 

C 

C  INPUT  NAME 

C  OUTPUT  ID, TYPE, DATE, Fl, SUB AV,T2,DATAB, NUMPT, Y 

C  NAME  =  CHR$  NAME  OF  THE  DATA  FILE 

C  ID    =  CHR$  TARGET  IDENTIFICATION 

C  TYPE  =  CHR$  WAVE  FORM  IDENTIFICATION 

C  DATE  =  CHR$  DATE 

C  INTEGER 

C  SUBAV  =  NUMBER  OF  SUB -AVERAGES  USED  TO  ACQUIRE  THE  DATA 

C  DATAB  =  DATA  BLOCK 

C  Fl    =  DUMMY 

C  NUMPT  =  NUMBER  OF  POINTS 

C  REAL 

C  T2    =  TIME  WINDOW 

C  Y     =  Y  ARRAY  OUTPUT 

CHARACTER''a6  NAME 
CHARACTER''-64  ID,F1 
CHARACTER''^64  TYPE, DATE 
REALM  Y(1500) 
INTEGER  SUBAV, DATAB, NUMPT, K 
PI=3. 141592654 
K=l 

OPEN  (UNIT=10,FILE=NAME,STATUS='OLD' ) 
C      READ  THE  HEADER 

READ(10,5)  ID 

READ(10,5)  TYPE 

READ(10,5)  DATE 

READ(10,5)  Fl 

READ( 10,2030)  SUBAV 

READ( 10,2010)  T2 

READ( 10,2030)  DATAB 

READ(10,20S0)  NUMPT 
10         READ(10,2000,END=20)  Y(K) 
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K=K+1 

17 

GOTO  10 

20 

CLOSE(IO) 

2000 

F0RMAT(E12.6) 

2010 

F0RMAT(F5. 2) 

2030 

F0RMAT(I5) 

5 

FORMAT(A) 

K=K-1 

IF(K. EQ. NUMPT)  GOTO  3000 

WRITE(''S*)  'CHECK  THE  DATA 

STOP 

3000 

RETURN 

END 

#  POINTS  MISMATCH' 


SUBROUTINE  FFF1(  NAME, YY,W,NEXPT, TIME) 


c 

WHERE 

* 

c 

DATE 

IS  THE  DATE  OF  MEASUREMENT 

* 

c 

DATAB 

IS  THE  NUMBER  OF  DATA  BLOCKS 

* 

c 

Fl 

IS  THE  DUMMY  STRING 

* 

c 

ID 

IS  THE  IDENTIFICATION  STRING 

* 

c 

K 

IS  THE  NUMBER  OF  POINT 

* 

c 

L9 

IS  THE  RELATIVE  SMOOTHING 

* 

c 

MODE 

IS  EQUAL  TO  0  MEAN  FORWARD 

■h 

c 

FOURIER  TRANSFORM 

c 

MAG 

IS  THE  MAGNITUDE  OF  Xo'VHt(f) 

* 

c 

NAME4 

IS  THE  NAME  OF  THE  TARGET  DATA 

iV 

c 

N 

IS  THE  NUMBER  OF  POINTS 

^V 

c 

NUMPT 

IS  THE  SAME  AS  K 

* 

c 

NEXPT 

IS  THE  POWER -OF -2  EXPONENT 

Vf 

c 

WILL  BE  GENERATED  BY  PREFFT 

•!V 

c 

SET  TO  -1  INDICATES  ERROR 

* 

c 

R 

IS  THE  REAL  PART  OF  THE  INVERSE 

* 

c 

FOURIER  TRANSFORM  OF  Xo'VHt(f) 

■k 

c 

SUBAV 

IS  THE  NUMBER  OF  SUB  AVERAGE 

•k 

c 

T 

IS  THE  TIME  INTERVAL 

* 

c 

TIME 

IS  THE  TIME  INTEVAL 

■>\ 

c 

TITLE 

IS  THE  TITLE  NAME  FOR  PLOITING 

Vf 

c 

TUT 

IS  THE  FREQUENCY  DOMAIN  OF  Xo*H(f)* 

c 

TYPE 

IS  THE  WAVEFORM  TYPE 

-'! 

c 

W 

IS  THE  COMPLEX  EXPONENTIAL  ARRAY 

>v 

c 

GENERATE  BY  PREFFT  SUBROUTINE 

•k 

c 

Xo*Ht(f) 

IS  THE  GAUSSIAN  WEIGHTED  COSINE 

it 

c 

XS 

IS  THE  STRING  FOR  PLOITING  ROUTINE* 

c 

FOR  THE  HEADER 

■k 

c 

DRIVER 

•k 

c 

Y4 

IS  THE  CALIBRATION  SPHERE 

■k 

c 
c 
c 

SUBTRACTED  WAVEFORM  (FREQUENCY) 

•k 

■k 

Y5 

IS  THE  TARGET  SUBTRACTED  WAVEFORM 

it 

c 
c 
c 

FREQUENCY  DOMAIN 

■k 
•k 

YSC 

IS  THE  COMPUTED  SCATTERED  ELECTRIC* 

c 

FIELD  FOR  THE  CANONICAL  TARGET 

it 

c 

Yl 

IS  THE  CALIBRATION  MEASUREMENT 

it 

c 

Y2 

IS  THE  BACKGROUND  MEASUREMENT 

it 
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c 

Y3 

IS  THE  TARGET  MEASURED 

* 

c 

YS 

IS  THE  SAME  AS  XS 

* 

c 

YN 

IS  THE  YES  OR  NO  STRING 

it 

c 

METER 

* 

c 

XMIl 

IS  THE  MINIMUM  VALUE 

* 

c 

XMAl 

IS  THE  MAXIMUM  VALUE 

* 

CHARACTER* 16  NAME 

CHARACTER''^64  ID, TYPE, DATE, TITLE, XS,YS,F1 

COMPLEX  W( 1025), YY( 1025) 

REAL  T 

REAL  Y( 1024) ,TIME ,T2 ,XMINS ,MAG( 1024) ,XMAG 

INTEGER  NUMPT,DATAB,SUBAV,MODE,N 

CALL  INPP(  NAME, ID, TYPE, DATE, F1,SUBAV,T2,DATAB,NUMPT,Y) 
TIME=T2/(FL0AT(NUMPT-l))'nE-9 
XMINS=0. 0 
N=NUMPT 

CALL  PRETIM( TITLE, XS,YS) 
CALL  P11(TITLE,N,XMINS,T2,Y,XS,ID,YS) 
C     SET  MODE  EQUA:  TO  0  TO  FIND  FOURIER  SERIES 
MODE=0 

CALL  PREFFT  (N,MODE,NEXPT,W) 
DO  1234  I  =1  ,NUMPT 

YY(I)=CMPLX(  Y(I),0. 0   ) 
1234         CONTINUE 
T=TIME 

CALL  FFT  (N,MODE,T,NEXPT,W,YY) 
C     CHANGING  THE  OUTPUT  YY  TO  MAGNITUDE  AND  PHASE 
N=NUMPT 

CALL  PREFRE( TITLE, XS,YS) 
CALL  MAGNITUD(YY,N,MAG) 
XMAG=1.  O/TIME/2.  O'nE-9 

CALL  P11(TITLE,N,XMINS,XMAG,MAG,XS,ID,YS) 
RETURN 
END 

Q*****************Vf  **********  Vr******yc****Vc****Vf************** 

C  * 

C  FFT  PROGRAM   FROM  THE  DIGITAL  SPECTRAL  ANALYSIS  * 

C  WITH  APLICATIONS   PRENTICE-HALL, INC. .  * 

C  ENGLEWOOD  CLIFF, NEW  JERSEY  07632  1987  * 

C  S.LAWRENCE  MARPLE  * 

C  * 

Q*********************************************************** 

SUBROUTINE  PREFFT  (N,MODE,NEXP,W) 
C 

C   Input  Parameters: 
C 

ON-  Number  of  data  samples  to  be  processed  (integer-must  be  a 
C  power  of  two) 

C     MODE  -  Set  to  0  for  discrete-time  Fourier  series  (Eq.  2.C. 1)  or 
C  1  for  inverse  (Eq.  2.  C.  2) 

C 

C     Output  Parameters: 
C 
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C     NEXT  -  Indicates  power-of-2  exponent  such  that  N=2*''''NEXP  . 

C  Will  be  set  to  -1  to  indicate  error  condition  if  N 

C  is  not  a  power  of  2  (this  integer  used  by  sub.  FFT) 

C     W     -  Complex  exponential  array 

C 

C     Notes: 

C 

C     External  array  W  must  be  dimensioned  .  GE.  N  by  calling  program. 

C 

COMPLEX  W(1024),C1,C2 

NEXP=1 
5     NT=2**NEXP 

IF  (NT  .GE.  N)   GO  TO  10 

NEXP=NEXP+1 

GO  TO  5 
10    IF  (NT  .EQ.  N)   GO  TO  15 

NEXP=-1 

RETURN 
15    S=8.  'VATAN(  1.  )/FLOAT(NT) 

C1=CMPLX(C0S(S),-SIN(S)) 

IF  (MODE  .NE.  0)   C1=C0NJG(C1) 

C2=(l.  ,0.) 

DO  20  K=1,NT 
W(K)=C2 
20      C2=C2*C1 

RETURN  - 

END 

SUBROUTINE  FFT  (N,MODE,T,NEXP,W,X) 
C 

0     Input  Parameters: 
C 

C     N,MODE,NEXP,V/  -  See  parameter  list  for  subroutine  PREFFT 
C     T  -  Sample  interval  in  seconds 

OX  -  Array  of  N  complex  data  samples,  X(l)  to  X(N) 

C 

C     Output  Parameters: 
C 

C     X  -  N  complex  transform  values  replace  original  data  samples 
C         indexed  from  k=l  to  k.=N,  representing  the  frequencies 
C        (k-l)/NT  hertz 
C 

C     Notes: 
C 

C     External  array  X  must  be  dimensioned  .  GE.  N  by  calling  program. 
C 

COMPLEX  X( 1024) ,W( 1024), CI, C2 

MM=1 

LL=N 

DO  70  K=1,NEXP 
NN=LL/2 
JJ=MM+1 

DO  40  I=1,N,LL 
KK=I+NN 
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C1=X(I)+X(KK) 
X(KK)=X(I)-X(KK) 
40        X(I)=C1 

IF  (NN  .EQ.  1)  GO  TO  70 
DO  60  J=2,NN 
C2=W(JJ) 
DO  50  I=J,N,LL 
KK=I+NN 
C1=X(I)+X(KK) 
X(KK)=(X(I)-X(KK))*C2 
50  X(I)=C1 

60       JJ=JJ+MM 
LL=NN 
MM=MM'>2 
70      CONTINUE 
NV2=N/2 
NM1=N-1 
J=l 
DO  90  1=1, NMl 

IF  (I  .GE.  J)   GO  TO  80 
C1=X(J) 
X(J)=X(I) 
X(I)=C1  . 
80      K=NV2 

85      IF  (K  .GE.  J)   GO  TO  90 
J=J-K 
K=K/2 
GO  TO  85 
90      J=J+K 

IF  (MODE  .EQ.  0)   S=T 
IF  (MODE  .NE.  0)   S=l. /(T^vflOAT(N)) 
DO  100  1=1, N 
100     X(I)=X(I)''-S 
RETURN 
END 

C 

SUBROUTINE  MAGNITUD  ( Y,NUMPT,MAG) 
C 

COMPLEX  Y(1024),H 
REAL*4  MAG(i024) 
INTEGER  N,K,NP,NUMPT 
DO  10  1=1  ,NUMPT/2 
H=Y(I) 

MAG(I)=CABS(H) 
10     CONTINUE 

NUMPT=NUMPT/2 

RETURN 

END 

SUBROUTINE  PRETIM( TITLE, XS.YS) 
CHARACTER'''64  TITLE ,  XS ,  YS 
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TITLE='TIME  SERIES  PLOT' 

XS   =' X-AXIS  TIME  IN  NS  ' 

YS   =' Y-AXIS  MAGNITUDE  IN  VOLT' 

RETURN 

END 

SUBROUTINE  PREFRE ( TITLE ,XS,YS) 

CHARACTER''^64  TITLE, XS,YS 

TITLE=' SPECTRAL  PLOT' 

XS=' X-AXIS  FREQUENCY  IN  GHZ' 

YS=' Y-AXIS  MAGNITUDE  IN  VOLT' 

RETURN 

END 

SUBROUTINE  IFFF(N,T,X) 

COMPLEX  W(1025),X(1025) 

REAL  T 

INTEGER  MODE 

M0DE=1 

CALL  PREFFT(N,MODE,NEXP,W) 

CALL  FFT(N,M0DE,T,NEXP,W3X) 

RETURN 

END 

SUBROUTINE  DEC0N3 (K, YSC, Y4, Y5 , ID, TYPE, DATE, Fl, SUB AV,T2,DATAB) 


c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 

COMPUTING  THE  TARGET  RESPONSE  TO  THE  INCIDENT 
WAVESHAPE  FOR  THE  COMPUTED  SPHERE.  USING 
RIAD'S  METHOD  OF  OPTIMAL  DECONVOLUTION. 

* 
* 

if 

Xo'VHt(f): 

=Y5*C0NJ(Y4)*Ysc/(Y4*C0NJ(Y4)+L) 

* 
■if 

WHERE 
Xo(f) 

IS  THE  INCIDENT  WAVEFORM 

■k 
* 
if 

Ht(f) 

IS  THE  TARGET  TRANSFER  FUNCTION 

if 

if 

Y4 

IS  THE  CALIBRATION  SPHERE 
SUBTRACTED  WAVEFORM  (FREQUENCY) 

* 
* 
if 

Y5 

IS  THE  TARGET  SUBTRACTED  WAVEFORM 
FREQUENCY  DOMAIN 

if 
* 
if 

YSC 

Yl 
Y2 
Y3 
K 

NUMPT 
R 

IS  THE  COMPUTED  SCATTERED  ELECTRIC 
FIELD  FOR  THE  CANONICAL  TARGET 
IS  THE  CALIBRATION  MEASUREMENT 
IS  THE  BACKGROUND  MEASUREMENT 
IS  TJ'E  TARGET  MEASURED 
IS  THE  NUMBER  OF  POINT 
IS  THE  SAME  AS  K 
IS  THE  REAL  PART  OF  THE  INVERSE 

* 
* 
* 
* 
it 
if 
if 
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C  FOURIER  TRANSFORM  OF  Xo*Ht(f)  * 

C  MAG  IS  THE  MAGNITUDE  OF  Xo*Ht(f)  * 

C  ID  IS  THE  IDENTIFICATION  STRING  * 

C  SUBAV  IS  THE  NUMBER  OF  SUB  AVERAGE  * 

C  DATAB  IS  THE  NUMBER  OF  DATA  BLOCK  * 

C  Fl  IS  THE  DUMMY  STRING  * 

C  NAME4  IS  THE  NAME  OF  THE  TARGET  DATA  * 

C  XS  IS  THE  STRING  FOR  PLOTTING  ROUTINE   * 

C  FOR  THE  HEADER  * 

C  YS  IS  THE  SAME  AS  XS  * 

C  TITLE  IS  THE  TITLE  NAME  FOR  PLOTTING  * 

C  YN  IS  THE  YES  OR  NO  STRING  * 

C  L9  IS  THE  SMOOTHING  PARAMETER  * 

C  TUT  IS  THE  FREQUENCY  DOMAIN  OF  Xo*H(f)    * 

C  TYPE  IS  THE  WAVEFORM  TYPE  * 

C  XMIl  IS  THE  MINIMUM  VALUE  * 

C  XMAl  IS  THE  MAXIMUM  VALUE  * 

C  DATE  IS  THE  DATE  OF  MEASUREMENT  * 


COMPLEX 

COMPLEX 

REAL 

REAL 

INTEGER 

CHARACTER* 16 

CHARACTER'>64 

CHARACTER'n 


Y5( 1025) ,Y4( 1025) , YSC( 1024) ,TUT( 1024) , SUM 

L1,L2,L3 

R( 1024), L9, XMAl, XMIl, MAG( 1024) 

TIME,T,T2 

K,NUMPT, SUBAV, DATAB 

NAME  4 

ID, YS,XS, TITLE, TYPE, DATE, Fl 

YN 

T=T2/FL0AT(K-l)*lE-9 

SUM=(0.  0,0.  0) 

DO  111  1=1, K 

SUM=SUM+Y4( I )*CONJG( Y4( I ) ) 
111   CONTINUE 
438   WR I TE  (  "  *  )  '  ■^■**'!f**->'f**if*-^'**-i-**-i-********-i'*-^******"***** ' 

WRITE(''V'0  '  ENTER  RELATIVE  SMOOTHING  PARAMETER  :     ' 

READ  (''-,*)  L9 

L1=CMPLX(L9,0. 0) 

L2=SUM/CMPLX(FL0AT(K+1) ,0. 0)*L1 

DO  1  1=1, K 

TUT(  I  )=Y5(  I) -'-CONJGC  Y4(I)  )*YSC(  I )/(  Y4(  I  )*CONJG(  Y4(  I )  )+L2) 
1     CONTINUE 

TITLE=' SPECTRAL  PLOT' 

YS='DECONVOLUTION' 

XS=' X-AXIS  FREQUENCY  IN  GHZ' 

NUMPT=K 

CALL  MAGNITUD(TUT,NUMPT,MAG) 

XMI1=0.  0 
XMAl=1.0/T/2.  0*lE-9 

CALL  P12(TITLE,NUMPT, XMIl, XMAl, MAG, XS, ID, YS,L9) 

CALL  IFFF(K,T,TUT) 

C'A-Vr***iViV'>ViV*')V*VrV(-VcVr*Vr*yr*****yr*****VrVf  **************** 

DO  222  1=1, K 
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R(I)=REAL(TUT(I)) 
222   CONTINUE 

XMI1=0.  0 

XMA1=T2 

XS=' X-AXIS  TIME  IN  NSEC' 

TITLE='TIME  SERIES  PLOT' 

YS='DECONVOLUTION' 

Q************************************************* 

CALL  P12(TITLE,K,XMI1,XMA1,R,XS,ID,YS,L9) 

C**********************'5V****Vr************Vc******** 

YN='N' 

WR I TE  ( *  *  )  *  ***********************•>>******************■*********** ' 

WRITE(*i''0  'WOULD  YOU  LIKE  TO  SAVE  THE  DECONVOLUTION  RESULT  ?:    ' 

READ(*,999)  YN 

IF  (  YN.  EQ.  'Y'.OR.  YN.  EQ.  'y')  GO  TO  4376 

GOTO  4388 

Q*************7V-!V******************************Vr*** 

4376  WRITEf'''^  ''') '  ****'''c*********Vc***********************' 

WRITE(''^',''0'  ENTER  FILE  NAME  FOR  THE  TIME  DOMAIN  DATA:   ' 

READ('V,999)  NAME4 

CALL  0UT(NAME4, ID, TYPE, DATE, Fl, SUB AV,T2, DAT AB,K,R) 
4388  YN='Y' 

WRITE  ( '"  '^^  ' **************************************************** 

WRITE  ('^*)  'REPEAT  WITH  NEW  SMOOTHING  PARAMETER  ?  (Y/N):  ' 

READ(*,999)  YN 

IF  (  YN.EQ.  'Y'.OR.  YN.EQ.  'y')  GO  TO  438 
999   FORMAT(A) 

RETURN 

END 

Q************************************************ 

C      F  =  ARRAY  OF  Y  VARIABLE 

C      XS=  $  X-SCALE  (STRING) 

C      YS=  $  Y-SCALE  (STRING) 

C      ID=  IDENTIFICATION  OF  THIS  GRAPH 

C 

SUBROUTINE  P  11( TITLE, NPTS,XMINN,X'MAXX,FFF,XS,  ID, YS) 
c 

C  RM-FORTRAN  Subroutine  to  Plot  a  Solid  Line  with  TITLE 

C  Based  on  Program  PLOT.  FOR.     Dec  1987  by  M.  A.  Morgan 

C 

C  INPUT  DATA  FORMAT 

C 

C  TITLE  -  64  Space  Title 

C  NPTS  -  #  Data  Points 

C  XMIN  -  Real  Min  X  value 

C  XMAX  -  Real  Max  X  value 

C  F(N)   -  Input  Data  Array 

C  IHC   -  Hardcopy  (l=Yes) 

C 

CHARACTER*!  YN,DUM 

CHARACTER*64  TITLE, XS, ID, YS 

REAL  X( 1025) ,F( 1025) , XMIN, XMAX, FFF( 1025) ,XMINN,XMAXX 

INTEGER*2  N , JROW , JCOL , I S YMBL , ITYPE , CYAN , WHITE , YELLOW 

IHC=0 
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WRITE (*,*)  '  WOULD  YOU  LIKE  TO  HAVE  A  HARD  COPY:  ' 

READ('^100)  YN 

IF(YN.  EQ.  'Y'.OR.  YN.  EQ.  'y')  IHC=1 

WHITE=7 
CYAN=11 
DO  1111  I=1,NPTS 
F(I)=FFF(I) 
1111   CONTINUE 

XMIN=XMINN 
XMAX=XMAXX 

YELL0W=14 

N=NPTS 

DX=(XMAX-XMIN)/(NPTS-1.  0) 

FMIN=0.  0 

FMAX=0.  0 

DO  22  K=1,NPTS 

X(K)=XMIN+(K-1.0)*DX 

IF(F(K).LT.  FMIN)  FMIN=F(K) 

IF(F(K).GT.  FMAX)  FMAX=F(K) 

22  CONTINUE 

IF(FMIN.  GT.  0.0)  FMIN=0.  0 
IF(FMAX.  LT.  0.0)  FMAX=0.  0 
!  Computing  Scale  Factors  for  Vertical  Axis 

ABSMIN=ABS(FMIN) 

ABSMAX=ABS(FMAX) 

YMAX=AMAX1(ABSMIN,ABSMAX) 

NSCL=INT(LOG10(YMAX)) 

IF  (YMAX.  LT.  1.  0)  NSCL=NSCL-1 

YSCL=10.**NSCL 

FMIN=FMIN/YSCL 

FMAX=FMAX/YSCL 

ABSMIN=ABSMIN/YSCL 

ABSMAX=ABSMAX/YSCL 

DO  33  K=1,NPTS 

33  F(K)=F(K)/YSCL 

YMIN=0. 0 
IF(FMIN.  EQ.  0.  0)  GO  TO  37 
35  YMIN=YMIN+0. 5 

IF(ABSMIN. GT. YMIN)  GO  TO  35 
YMIN=YMIN''-FMIN/ABSMIN 
37  CONTINUE 

YMAX=0. 0 
IF(FMAX.  EQ.  0.0)  GO  TO  41 
39  YMAX=YMAX+0. 5 

IF(ABSMAX.  GT. YMAX)  GO  TO  39 
YMAX=YMAX*FMAX/ AB SMAX 

41  CONTINUE 

CALL  QBEEP 

CALL  QSM0DE(2) 
WRITE(*,*)  'Press  RET  to  View  Plot.  ..  RET  Again  to  Clear  Screen' 

IF(IHC.NE.  1)  GO  TO  42 
WRITE(*,*)  'A  Hardcopy  Will  Be  Made. .. Check  that  Printer  is  ON' 

42  CONTINUE 

READ(*,100)  DUM 

:  Calling  GRAFMATIC  Routines  and  Plotting  Solid  Line  Graph 

ITYPE=1 


57 


ISYMBL=-2 
CALL  QSM0DE(16) 
CALL  QPTXTC 64, TITLE, YELLOW, 16,23) 
CALL  QPTXT(64,ID,CYAN,16,19) 
CALL  QPTXTC 64, XS, CYAN, 16,20) 
CALL  QPTXTC 64, YS,CYAN, 16,21) 

CALL  QPLOTC110,540,20,240,XMIN,XMAX,YMIN,YMAX,0.  ,0.  ,0,1.  ,1.5) 
CALL  QSETUPCO,CYAN,ISYMBL,CYAN) 
XMAJ0R=XMAX/5.0 

CALL  QXAXISCXMIN,XMAX,XMAJOR, 1,1,2) 
YMAJ0R=1.  0 

CALL  QYAXISCYMIN,YMAX,YMAJOR, 1,1,2) 

JROW=14+CABSCYMIN)/CABSCYMAX)+ABSCYMIN)))*220 

JCOL=80-430*XMIN/CXMAX-XMIN) 

CALL  QGTXTC3, '0.  0' , WHITE, JCOL,JROW,0) 

CALL  QPTXTCl,'s' ,YELL0W,5,16) 

CALL  QPTXTC 1 , ' c ' , YELLOW ,5,15) 

CALL  QPTXTCl,'a' ,YELL0W,5,14) 

CALL  QPTXTCl,'!' ,YELL0W,5,13) 

CALL  QPTXTC 1, ' e' , YELLOW, 5, 12) 

CALL  QPTXTC 1 , ' X '  YELLOW ,5,10) 

CALL  QPTXTC 2,' 10 ', YELLOW, 4, 8) 

CALL  QPTXTC 1 , ' * ' , YELLOW ,5,7) 

CALL  QPTXTCl,'*' ,YELL0W,5, 6) 

CALL  QCMOVC0,6) 

WRITEC*,150)  NSCL 

CALL  QTABLCITYPE,N,X,F) 

READC*,100)  DUM 

IFCIHC.NE.  1)  GO  TO  44 

CALL  PRTSC 

44  CONTINUE 

CALL  QSM0DEC2) 
100  FORMATCA) 

150  F0RMATC4X,I3) 

RETURN 
END 

0 

C  THIS  SUBROUTINE  PLOTS  TWO  GRAPHS  ON  THE  SAME  AXIS 

C 

C  F  =  Yl  ARRAY  VARIABLE  C  BLUE) 

C  Fl  =  Y2  ARRAY  VARIABLE  C GREEN) 

C  AFTER  PLOT  THIS  PROGRAM  BE  AWARE  OF  THE  F1,F 

C  IT  WILL  SENT  BACK  A  SCALE  F,F1 

C 

C  XS  =  $  X- SCALE 

C  YS  =  $  Y-SCALE 

C  ID  =  IDENTIFICATION 

SUBROUTINE  PL08CTITLE,NPTS,XMINN,XMAXX,PUM,XS,YS,PUM1,ID) 
C 

C  RM-FORTRAN  Subroutine  to  Plot  a  Solid  Line  with  TITLE 

C  Based  on  Program  PLOT. FOR.    Dec  1987  by  M.  A.  Morgan 

C 
C  INPUT  DATA  FORMAT 
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c 

C  TITLE  -  64  Space  Title 

C  NPTS   -  #  Data  Points 

C  XMIN   -  Real  Min  X  value 

C  XMAX   -  Real  Max  X  value 

C  F(N)   -  Input  Data  Array 

C  IHC   -  Hardcopy  (l=Yes) 

C 

CHARACTER*!  YN,DUM 

CHARACTER'>64  TITLE, XS, ID, YS 

REAL  X( 1025) ,F( 1025) ,F1( 1025) ,PUM( 1025) ,PUM1( 1025) ,XMINN,XMAXX 

INTEGER*2  N , JROW , JCOL , I S YMBL , ITYPE , CYAN .WHITE , YELLOW , GREEN 

WRITE (*,*)  'WOULD  YOU  LIKE  TO  HAVE  A  HARD  COPY  ?:   ' 
READ(*,100)  YN 
IHC=0 

IF(YN.  EQ.  'Y'.OR.  YN.  EQ.  'y')  IHC=1 
DO  159  1=1, NPTS 
F(I)=PUM(I) 
F1(I)=PUM1(I) 
159    CONTINUE 

XMIN=XMINN 
XMAX=XMAXX 

WHITE=7 
CYAN=11 
YELL0W=14 
GREEN=2 

N=NPTS 

DX=(XMAX-XMIN)/(NPTS-1.  0) 

FMIN=0.  0 

FMAX=0. 0 

FMIN1=0.  0 

FMAX1=0.  0 

DO  22  K=1,NPTS 
X(K)=XMIN+(K-1.0)*DX 
IF(F(K).LT.  FMIN)  FMIN=F(K) 
IF(F(K).GT. FMAX)  FMAX=F(K) 
IF(F1(K).LT.  FMINl)  FMIN1=F1(K) 
IF(F1(K).GT.  FMAXl)  FMAX1=F1(K) 
22    CONTINUE 

FMIN=AMIN1( FMIN, FMINl) 
FMAX=AMAX 1 ( FMAX , FMAX 1 ) 

IF(FMIN.  GT.  0.  0)  FMIN=0.  0 
IF(FMAX.  LT.  0.0)  FMAX=0.  0 
C  Computing  Scale  Factors  for  Vertical  Axis 

ABSMIN=ABS(FMIN) 

ABSMAX=ABS(FMAX) 

YMAX=AMAX1( ABSMIN , ABSMAX) 

NSCL=INT(LOG10(YMAX)) 

IF  (YMAX.  LT.  1.0)  NSCL=NSCL-1 

YSCL=10.**NSCL 

FMIN=FMIN/YSCL 

FMAX=FMAX/YSCL 

ABSMIN=ABSMIN/YSCL 

ABSMAX=ABSMAX/YSCL 
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DO  33  K=1,NPTS 
F1(K)=F1(K)/YSCL 
33  F(K)=F(K)/YSCL 

YMIN=0. 0 
IF(FMIN.  EQ.  0.0)  GO  TO  37 
35  YMIN=YMIN+0.  5 

IFCABSMIN. GT. YMIN)  GO  TO  35 
YMIN=YMIN'VFMIN/ABSMIN 
37  CONTINUE 

YMAX=0. 0 
IF(FMAX.  EQ.  0.0)  GO  TO  41 
39  YMAX=YMAX+0. 5 

IF(ABSMAX. GT.  YMAX)  GO  TO  39 
YMAX= YMAX*FM AX / AB  SMAX 

41  CONTINUE 

CALL  QBEEP 

CALL  QSM0DE(2) 
WRITE(''S*)  'Press  RET  to  View  Plot... RET  Again  to  Clear  Screen' 

IF(IHC.NE.  1)  GO  TO  42 
WRITE(*,*)  'A  Hardcopy  Will  Be  Made. .. Check  that  Printer  is  ON' 

42  CONTINUE 

READC'slOO)  DUM 

C  Calling  GRAFMATIC  Routines  and  Plotting  Solid  Line  Graph 

ITYPE=1 
ISYMBL=-2 
CALL  QSM0DE(16) 
CALL  QPTXTC  64 , TITLE , YELLOW ,16,23) 
CALL  QPTXT(64,ID,GREEN,16,19) 
CALL  QPTXT(64,XS, WHITE, 16,20) 
CALL  QPTXTC 64, YS, CYAN, 16,21) 

CALL  QPLOTC 110, 540, 20, 240, XMIN,XMAX, YMIN, YMAX, 0.  ,0.  ,0,1.  ,1.5) 
CALL  QSETUPCO, CYAN, ISYMBL, CYAN) 
XMAJOR=XMAX/5.0 
CALL  QXAXISfXMIN,XMAX,XMAJOR, 1,1,2) 

YMAJ0R=1. 0 

CALL  QYAXIS(YMIN,YMAX,YMAJOR, 1,1,2) 

JR0W=14+(ABS( YMIN) /(ABS(YMAX)+ABS( YMIN) ))*220 

JC0L=80-430''-XMIN/(XMAX-XMIN) 

CALL  QGTXT(3,'0.0' , WHITE, JCOL,JROW,0) 

CALL  QPTXT(l,'s' ,YELLOW,5,16) 

CALL  QPTXTC 1, ' c' , YELLOW, 5, 15) 

CALL  QPTXTC 1,' a' , YELLOW, 5, 14) 

CALL  QPTXTC l,'l', YELLOW, 5, 13) 

CALL  QPTXTC 1, ' e' , YELLOW, 5, 12) 

CALL  QPTXTC 1 , ' X ' , YELLOW ,5,10) 

CALL  QPTXTC2,'10' ,YELL0W,4,8) 

CALL  QPTXTC 1 , ' * ' , YELLOW ,5,7) 

CALL  QPTXTC 1, '*' ,YELL0W,5, 6) 

CALL  QCM0VC0,6) 

WRITEC'-,150)  NSCL 

CALL  QTABLCITYPE,N,X,F) 

C     PLOT  ANOTHER  GRAPH 

GREEN=10 

ISYMBL=-2 
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ITYPE=1 

CALL  QSETUP(0, GREEN, I SYMBL, GREEN) 

CALL  QTABL(ITYPE,N,X,F1) 

READ('^100)  DUM 

IF(IHC.NE.  1)  GO  TO  44 

CALL  PRTSC 

44  CONTINUE 

CALL  QSM0DE(2) 
100  FORMAT(A) 

150  F0RMAT(4X,I3) 

RETURN 
END 

C  F  =  ARRAY  OF  Y  VARIABLE 

C      XS=  $  X-SCALE  (STRING) 

C      YS=  $  Y-SCALE  (STRING) 

C      ID=  IDENTIFICATION  OF  THIS  GRAPH 

C 

SUBROUTINE  P 12( TITLE, NPTS,XMINN,XMAXX,FFF,XS, ID, YS,L9) 
c 

C  RM-FORTRAN  Subroutine  to  Plot  a  Solid  Line  with  TITLE 

C  Based  on  Program  PLOT. FOR.     Dec  1987  by  M.  A.  Morgan 

C 

C  INPUT  DATA  FORMAT 

C 

C  TITLE  -  64  Space  Title 

C  NPTS   -  #  Data  Points 

C  XMIN  -  Real  Min  X  value 

C  XMAX  -  Real  Max  X  value 

C  F(N)   -  Input  Data  Array 

C  IHC    -  Hardcopy  (l=Yes) 

C 

CHARACTER* 1  YN,DUM 

CHARACTER*64  TITLE, XS, ID, YS ,RELAT 

REAL  X( 1025) ,F( 1025) , XMIN, XMAX, FFF( 1025) ,XMINN,XMAXX,L9 

INTEGER*2  N , JROW , JCOL , I SYMBL , ITYPE , CYAN , WHITE , YELLOW 

IHC=0 

WRITE(*,*)  '  WOULD  YOU  LIKE  TO  HAVE  A  HARD  COPY:  ' 

READ(*,100)  YN 

IF(YN.  EQ.  'Y'.OR.  YN.  EQ.  'y')  IHC=1 

WHITE=7 

CYAN=11 

DO  1111  1=1, NPTS 


F(I)=FFF(T) 
1111   CONTINUE 


XMIN=XMINN 

XMAX=XMAXX 

YELL0W=14 

N=NPTS 

DX=(XMAX-XMIN)/(NPTS-1.  0) 

FMIN=0. 0 

FMAX=0.  0 

DO  22  K=1,NPTS 

X(K)=XMIN+(K-1.0)*DX 

IF(F(K).LT.  FMIN)  FMIN=F(K) 
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IF(F(K).GT.  FMAX)  FMAX=F(K) 
22  CONTINUE 

IF(FMIN.  GT.  0.  0)  FMIN=0.  0 

IF(FMAX.  LT.  0.  0)  FMAX=0.  0 

Computing  Scale  Factors  for  Vertical  Axis 

ABSMIN=ABS(FMIN) 

ABSMAX=ABS(FMAX) 

YMAX=AMAX1(ABSMIN,ABSMAX) 

NSCL=INT(LOG10(YMAX)) 

IF  (YMAX.  LT.  1.0)  NSCL=NSCL-1 

YSCL=10.  **NSCL 

FMIN=FMIN/YSCL 

FMAX=FMAX/YSCL 

ABSMIN=ABSMIN/YSCL 

ABSMAX=ABSMAX/YSCL 

DO  33  K=1,NPTS 

33  F(K)=F(K)/YSCL 

YMIN=0.  0 
IF(FMIN.  EQ.  0.  0)  GO  TO  37 
35  YMIN=YMIN+0.  5 

IF(ABSMIN. GT. YMIN)  GO  TO  35 
YMIN=YMIN*FMIN/ABSMIN 
37  CONTINUE 

YMAX=0.  0 
IF(FMAX.  EQ.  0.  0)  GO  TO  41 
39  YMAX=YMAX+0. 5 

IF(ABSMAX. GT. YMAX)  GO  TO  39 
YMAX=YMAX'''FMAX/ABSMAX 

41  CONTINUE 

CALL  QBEEP 

CALL  QSM0DE(2) 
WRITE(*,*)  'Press  RET  to  View  Plot.  ..RET  Again  to  Clear  Screen' 

IF(IHC.NE.  1)  GO  TO  42 
WRITEC'S")  'A  Hardcopy  Will  Be  Made. .. Check  that  Printer  is  ON' 

42  CONTINUE 

READ(''-,100)  DUM 

Calling  GRAFMATIC  Routines  and  Plotting  Solid  Line  Graph 

ITYPE=1 

ISYMBL=-2 

CALL  QSM0DE(16) 

CALL  QPTXT( 64 , TITLE , YELLOW ,16,23) 

CALL  QPTXT(64,ID,CYAN,16,19) 

CALL  QPTXT( 64, XS, CYAN, 16,20) 

CALL  QPTXT( 64, YS, CYAN, 16,21) 

RELAT=' RELATIVE  SMOOTH  =' 

CALL  QPTXT(64,RELAT,CYAN,16,18) 

CALL  QQNPUT(265,253,L9,3) 

CALL  QPLOT( 110, 540, 20, 240, XMIN,XMAX, YMIN, YMAX, 0.  ,0.  ,0,1.  ,1.5) 

CALL  QSETUP( 0, CYAN, ISYMBL, CYAN) 

XMAJ0R=XMAX/5.  0 

CALL  QXAXIS(XMIN,XMAX,XMAJOR, 1,1,2) 

YMAJ0R=1.  0 

CALL  QYAXIS(YMIN,YMAX,YMAJOR, 1,1,2) 

JROW=14+(ABS(YMIN)/(ABS(YMAX)+ABS(YMIN)))*220 

JC0L=80-430'''XMIN/(XMAX-XMIN) 

CALL  QGTXTC  3 , ' 0 . 0 ' , WHITE , JCOL , JROW , 0 ) 
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CALL  QPTXT(l,'s' ,YELLOW,5,16) 

CALL  QPTXT( 1 , ' c ' , YELLOW ,5,15) 

CALL  QPTXT( 1 , ' a ' , YELLOW ,5,14) 

CALL  QPTXT(l,'r  ,YELL0W,5,13) 

CALL  QPTXT(l,'e' ,YELL0W,5,12) 

CALL  QPTXT( 1 , ' X '  YELLOW ,5,10) 

CALL  QPTXT(2,'10' ,YELL0W,4,8) 

CALL  QPTXT(1,'*',YELL0W,5,7) 

CALL  QPTXTd,'*' ,YELL0W,5,6) 

CALL  QCM0V(0,6) 

WRITE(*,150)  NSCL 

CALL  QTABL(ITYPE,N,X,F) 

READ(*,100)  DUM 

IFCIHC.NE.  1)  GO  TO  44 

CALL  PRTSC 

44  CONTINUE 

CALL  QSM0DE(2) 
100  FORMAT(A) 

150  F0RMAT(4X,I3) 

RETURN 
END 

SUBROUTINE  PL077(Y1 ,Y5,T2,NUMPT, IDl) 

C        Y2  =  Y5  =  BACKGROUND  DATA  ARRAY  * 

C        Yl      =  SCATTERED  DATA  ARRAY  * 

INTEGER  NUMPT,DATAB,SUBAV,F1,NSHIF,SET 

REAL  Yl( 1024) , Y2( 1024) , Y3( 1024) , Y4( 1024) ,Y5( 1024) 

REAL  XMIN,XMAX 

CHARACTER^ne  NAME  1 , NAME 2, NAME 3, QN 

CHARACTER'n  YN 

CHARACTER''^64  YS,XS,  ID, TITLE, DATE 

CHARACTER''^64  ID1,TYPE1 

CHARACTER'^64  ID2,TYPE2 

CHARACTER'''64  ID3,TYPE3 

REAL  TIME,T2 

WR I TE  ( '"  " )  '  * Vf Vf Vf Vf V-** ff Vf vc******Vf*yf ****y-*yf ?v vciv*****Vf*********Vc**** A* ' 

WRITE('^'V)  'SHIFTING  SCATTERED  WAVEFORM  &  SUBTRACTING  BACKGROUND' 

XMIN=0. 6 
XMAX=T2 

DO  22  K=1,NUMPT 
22  Y2(K)=Y5(K) 
SET=0 

Q****VfV-yfyr***VfycVf**VfVf****y('**':lf*yf***yf**Vf  ************* 

33  CALL  SHIFT(TITLE,YS,XS,ID) 

Q***Vr************yf  ******************************* 

CALL  PL08(TITLE,NUMPT,XMIN,XMAX,Y1,XS,YS,Y2,ID) 
C**********************Vf****Vf******************** 

C     SHIFT  THE  CURVE 
C************************************************ 

CALL  NSHIFT(NSHIF) 
IFCNSHIF.EQ.  0)  GO  TO  40 
CALL  G00D(Y1,Y2,NSHIF,NUMPT,SET) 
GO  TO  33 


63 


40  CONTINUE 
C     SUBTRACT  BACKGROUND  FROM  THE  WAVEFORM  **** 

DO  44  I=1,NUMPT 
44  Y3(I)=Y1(I)-Y2(I) 

CALL  SUBTR( TITLE, XS,YS) 

CALL  P11(TITLE,NUMPT,XMIN,XMAX,Y3,XS,ID1,YS) 

WRITE(*,'V)  'DO  YOU  WANT  TO  RE-SHIFT  AND  SUBTRACT  ?  (Y/N): 

READ(''-,300)  YN 

IF(YN.  EQ.  'y'.OR.  YN.  EQ.  'Y')  GO  TO  33 

DO  77  I=1,NUMPT 
77  Y1(I)=Y3(I) 
300  FORMAT(A) 

RETURN 

END 

SUBROUTINE  NSHIFT(NSHIF) 
INTEGER  NSHIF 

WR I TE  (  '"'■  ''^  )   '  '"f  ^'V'jViVVriV Vf ***VciV*****iV*':^*'jV******Vryr*Vr**'»VVf *****  ' 

WRITE(''S''0 

WRITE(''S'V)  'SCATTERED  WAVEFORM  SHIFTING:   ' 

WRITE  (''%-''') 

WRITE(",'''0  '+n  shifts  scat  waveform  to  the  RIGHT  n-points' 

WRITER'S*)  ' -n  shifts  scat  waveform  to  the  LEFT  n-points' 

WRITE(''S''0  '  0  DISPLAYS  subtracted  waveform   ' 

WRITE  (*,''0 

WRITE('"S''0  'ENTER  SHIFT  (+/-INTEGER):   ' 

READ(-'V,''0  NSHIF 

RETURN 

END 

C***************************************Vf  ******** 

SUBROUTINE  SUBTR( TITLE, XS,YS) 

CHARACTER-'64  TITLE, XS,  ID, YS 

TITLE=' SUBTRACTED  WAVEFORM' 

XS=' X-AXIS  TIME  IN  NS' 

YS=' Y-AXIS  MAGNITUDE  IN  VOLT' 

RETURN 

END 

(]*****•****************************•*************** 

SUBROUTINE  SHIFT(TITLE,YS,XS,ID) 

CHARACTER*64  TITLE, YS,XS, ID 

TITLE='TIME  SERIES  PLOT' 

YS=' BLUE -TARGET' 

ID= ' GREEN -BACKGROUND ' 

XS=' SHIFT-TARGET  ONLY' 

RETURN 

END 

(]*****************Vf*i'.-5V******v-****************yf**** 
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SUBROUTINE  G00D(Y1,Y2,NSHIF,NUMPT,SET) 

C  * 

C  PROGRAM  TO  SHIFT  Y2  CURVE  * 
C 

C  -   SHIFT  Yl  TO  THE  LEFT  * 

C  +   SHIFT  Yl  TO  THE  RIGHT  * 
C 

C  Yl       IS  THE  SHIFTING  CURVE  * 

C  Y2       IS  THE  BACKGROUND  CURVE  * 

C  NSHIF    IS  THE  NUMBER  OF  POINT  TO  SHIFT  * 

C  SET      IS  THE  SUM  OF  THE  POINT  TO  SHIFT  * 

C  IN  ORDER  TO  ZERO  PADDING  * 

C  NUMPT    IS  THE  NUMBER  OF  POINT  * 

C  * 

REAL  Yl( 1024), Y2( 1024), Y3(1024) 
INTEGER   NSHIF, NUMPT, SET 
5    IF(  NSHIF.  LT.  0  )  THEN 

NSHIF=-1*( NSHIF) 
DO  20  I=NSHIFT,NUMPT-NSHIF-SET 
Y3(I)=Y1(I+NSHIF) 
20  CONTINUE 

DO  10  I=NUMPT-NSHIF- SET, NUMPT 
Y2(I)=0.0 
Y3(I)=0.0 
10  CONTINUE 

DO  15  I=1,NSHIF+SET+1 
Y2(I)=0.0 
Y3(I)=0.  0 
15  CONTINUE 

ELSE  IF(  NSHIF. GT. 0  )  THEN 
N=NSHIF+SET+1 
DO  40  I=N, NUMPT 

Y3(I)=Y1(I-NSHIF) 
40  CONTINUE 

DO  30  1=1, N 

Y2(I)=0.0 
Y3(I)=0. 0 
30  CONTINUE 

DO  35  I=NUMPT-NSHIF-SET, NUMPT 
Y2(I)=0.0 
Y3(I)=0.0 
35  CONTINUE 

ELSE 
END  IF 

SET=SET+NSHIF 
DO  50  1=1, NUMPT 
Y1(I)=Y3(I) 
50    CONTINUE 
RETURN 
END 

SUBROUTINE  OUT(NAME , ID,TYPE , DATE ,F1 ,SUBAV,T2 ,DATAB , NUMPT, Y) 
C 
C     ID  $     =  CHARACTER  $  TARGET  IDENTIFICATION 
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c 

TYPE  $ 

c 

DATE  $ 

c 

Fl 

c 

SUBAV 

c 

NUMPT 

c 

DATAB 

c 

Y 

G 

T2 

=  TYPE  OF  WAVEFORM 

=  DATE  OF  EXPERIMENT  TAKE  PLACE 

=  DUMMY  STRING 

=  SUB  AVERAGE 

=  NUMPT  OF  POINT 

=  DATA  BLOCK 

=  Y  ARRAY 

=  TIME  WINDOW 
CHARACTER'>64  ID, TYPE, DATE, Fl 
CHARACTER'ne  NAME 
REAL  Y(1024),T2 
INTEGER  SUBAV, NUMPT, DATAB 
OPEN( UNIT=10 ,FILE=NAME , STATUS=' UNKNOWN' ) 
WRITE(10,'V)  ID 

WRITEdO,*)  TYPE 
WRITECIO,*)  DATE 
WRITEdO,*)  Fl 
WRITE( 10,22)  SUBAV 
WRITE( 10,21)  T2 
WRITE( 10,22)  DATAB 
WRITE( 10,22)  NUMPT 

DO  10  1=1, NUMPT 

WRITE( 10,20)  Y(I) 
10  CONTINUE 

20  F0RMAT(E12.6) 

21  FORMAT(F5.2) 

22  F0RMAT(I5) 
30   FORMAT(A) 

CLOSE(IO) 

RETURN 

END 
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APPENDIX  B.      THE  ACQUISITION  QUICK  BASIC  PROGRAM 

The  following  computer  program  acquires  the  measurement  data  from  the  HP 
54120T  as  described  in  chapter  5.   The  program  is  written  in  QUICK  BASIC 

1000  REM  ********************* ********Vc***Vc**Vf****Vr* 

1020  REM  *   MODIFY  4/22/88  * 

1040  REM  *       WRITTEN  BY  LT. JG  SOONPUEN  SOMAPEE  * 

1060  REM  *   ACQUISITION  * 

1080  REM  *  * 

1240  REM  *  * 

1321  REM  ******************************************* 

1322  D$=DATE$ 

1324  DUMMY$=TIME$:  PRINT  "DATE ";  D$:  PRINT  "TIME ";  DUMMY$ 

1325  COLOR  7,1 

1326  PRINT  "WELCOME  TO  ACQUISITION  PROGRAM" 

1327  PRINT  "THERE  ARE  3  TYPE  OF  MEASUREMENT" 

1328  PRINT  "1.  BACKGROUND WAVETYPE" 

1329  PRINT  "2.  CALIBRATE WAVETYPE" 

1330  PRINT  "3.  TARGET WAVETYPE" 

1341  PRINT  "  PRESS  THE  FOLLOWING  KEY  " 

1342  PRINT  "  L    (  USED  DATA  LAST  RESPONSE) " 

1343  PRINT  "  D    (  SYSTEM  DEFAULT  ) " 

1344  V$=INKEY$: IF  V$=""  THEN  1344 

1345  IF  V$="L"  OR  V$  ="l"  THEN  1960 

1346  IF  V$="d"  OR  V$  ="D"  THEN  1360 

1347  GOTO  1344 

1360  REM  ******************************************* 
1380  REM  ^^    DEFAULT  ^' 

1400  REM  ******************************************* 

1420  OPEN  "DATA. DAT"  FOR  OUTPUT  AS  #1 

1440   Pl=512 

1460   PRINT  #1,P1 

1480   N=1024 

1500   PRINT#1,N 

1520   ND=1 

1540   PRINT// 1,ND 

1560  T$="4   inch  Thin  wire" 

1580   PRINT// 1,T$ 

1600  D$=DATE$ 

1620  PRINT//1,D$ 

1640  T2=1.563E-10 

1660  PRINT//1,T2 

1680  W$="Background  Wave  Type" 

1700  PRINT/ZljW? 

1720  Fl=20 

1740  PRINT#1,F1 

1760  Sl=16 

1780  PRINT//1,S1 

1800  G1=S1 

1820  DSE$="Backg" 

1840  PRINT// 1,DSE$ 
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1860 
1880 
1900 
1920 
1940 
1960 
1980 
2000 
2020 
2040 
2060 
2080 
2100 
2120 
2140 
2160 
2180 
2200 
2220 
2240 
2260 
2280 
2300 
2320 
2340 
2360 
2380 
2400 
2420 
2440 
2460 
2480 
2500 
2520 
2540 
2560 
2580 
2600 
2620 
2640 
2660 
2680 
2700 
2720 
2740 
2760 
2780 
2800 
2820 
2840 
2860 
2880 
2900 
2920 
2940 
2960 


DUMMY$=TIME$ 

PRINT/n,DUMMY$ 

DELAY1=38 

PR  I  NT// 1,  DELAY  1 

CLOSED 1 

REM  Optionally  include  the  following  declarations 

REM  in  your  program. 

REM  They  provide  appropriate  mnemonics  by  which 

REM  to  reference  commonly  used  values. 

REM  Some  mnemonics  (GET%,  ERR%, 

REM  END%,  ATN%)  are  preceded  by  "B"  in  order  to 

REM  distinguish  them  from 

REM  BASICA  keywords. 

REM 


REM  GPIB  Commands 
UNL%  =  &H3F 
UNr^  =  &H5F 
GTL%  =  &H1 
SDC%  ^  &H4 
PPC%  =  &H5 
BGET%  =  &H8 
TCT%  =  &H9 
LLO%  =  &H11 
DCL%  =  6cH14 
PPU%  =  6cH15 
SPE%  =  &H18 
SPD%  =  &H19 
PPE%  =  ScH60 
PPD%  =  5cH70 


GPIB  unlisten  command 
GPIB  untalk  command 
GPIB  go  to  local 
GPIB  selected  device  clear 
GPIB  parallel  poll  configure 
GPIB  group  execute  trigger 
GPIB  take  control 
GPIB  local  lock  out 
GPIB  device  clear 
GPIB  ppoll  unconfigure 
GPIB  serial  poll  enable 
GPIB  serial  poll  disable 
GPIB  parallel  poll  enable 
GPIB  parallel  poll  disable 


REM 

REM  GPIB  status  bit  vector 

REM  global  variable  IBSTA%  and  wait  mask 

Error  detected 
Timeout 

EOI  or  EOS  detected 
SRQ  detected  by  CIC 
Device  needs  service 
I/O  completed 
Local  lockout  state 
Remote  state 
Control ler-In-Charge 
Attention  asserted 
Talker  active 
Listener  active 
Device  trigger  state 
Device  clear  state 


BERR% 

=  &H8000   ' 

TIMOX 

=  &H4000   ' 

BEND% 

=  &H2000   ' 

SRQI% 

^  &H1000   ' 

RQS%  = 

=  &Hc:00 

CMPL% 

=  &H100    ' 

LOK%  = 

=  &H80      • 

REM%  = 

:  &H40      ' 

CIC%  = 

=  &H20 

BATN% 

=  &H10    ' 

TACS% 

=  &H8 

LACS% 

=  &H4     • 

DTAS% 

=  &H2     • 

DCAS% 

=  &H1     • 

REM 

REM  Error  message 

EDVR% 

=  0 

ECIC% 

^   1      • 

EN0L7o 

=  2      ' 

EADR% 

=  3      ' 

EARG% 

==  4      ' 

ESAC% 

=  5 

EABO% 

=  6 

DOS  error 

Function  requires  GPIB-PC  to  be  CIC 
Write  function  detected  no  Listeners 
Interface  board  not  addressed  correctly 
Invalid  argument  to  function  call 
Function  requires  GPIB-PC  to  be  SAC 
I/O  operation  aborted 
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2980 
3000 
3020 
3040 
3060 
3080 
3100 
3120 
3140 
3160 
3180 
3200 
3220 
3240 
3260 
3280 
3300 
3320 
3340 
3360 
3380 
3400 
3420 
3440 
3460 
3480 
3500 
3520 
3540 
3560 
3580 
3600 
3620 
3640 
3660 
3680 
3700 
3720 
3740 
3760 
3780 
3800 
3820 
3840 
3860 
3880 
3900 
3920 
3940 
3960 
3980 
4000 
4020 
4040 
4060 
4080 


ENEB%  = 

EOIP%  = 

REM 

ECAP%  = 

EFSO%  = 

EBUS%  = 

ESTB%  = 

ESRQ%  = 

REM 

REM  EOS 

BIN%  =  ScHlOOO 


7 
10 

11 
12 
14 
15 
16 


Non-existent  interface  board 

I/O  operation  started  before  previous 

operation  completed 

No  capability  for  operation 

File  system  operation  error 

Command  error  during  device  call 

Serial  poll  status  byte  lost 

SRQ  remains  asserted 


mode  bits 


XEOS%  =  &H800 

REOS%  =  &H400 

REM 

REM  Timeout  values 


Eight  bit  compare 
Send  EOI  with  EOS  byte 
Terminate  read  on  EOS 


TNONE%  = 
T10US%  = 
T30US%  = 
T100US%  = 
T300US%  = 


0 
1 
2 
3 

4 
5 
6 

7 
8 
9 
10 
11 
12 
13 
14 
15 
16 
17 


and  meanings 
Infinite  timeout  (disabled) 
Timeout  of  10  us  (ideal) 
Timeout  of  30  us  (ideal) 
Timeout  of  100  us  (ideal) 
Timeout  of  300  us  (ideal) 
Timeout  of  1  ms  (ideal) 
Timeout  of  3  ras  (ideal) 
Timeout  of  10  ms  (ideal) 
Timeout  of  30  ms  (ideal) 
Timeout  of  100  ms  (ideal) 
Timeout  of  300  ms  (ideal) 
Timeout  of  1  s  (ideal) 
Timeout  of  3  s  (ideal) 
Timeout  of  10  s  (ideal) 
Timeout  of  30  s  (ideal) 
Timeout  of  100  s  (ideal) 
Timeout  of  300  s  (ideal) 
Timeout  of  1000  s  (maximum) 


T1MS%  = 

T3MS%  = 

T10MS%  = 

T30MS%  = 

T100MS%  = 

T300MS%  = 

T1S%  = 

T3S%  = 

T10S%  = 

T30S%  = 

T100S%  = 

T300S%  = 

T1000S%  = 

REM 

REM  Miscellaneous 

S%  =  &H8 

LF%  =  &HA       • 

REM 

REM  Application  program  variables  passed  to 

REM  GPIB  functions 

REM 

CMD$  =  SPACE$(10) 

RD$  =  SPACE$(255) 

WRT$  =  SPACE$(255) 

BNAME$  =  SPACE$(7) 

BDNAME$  =  SPACE$(7) 

FLNAME$  =  SPACE$(50) 


Parallel  Poll  sense  bit 
Line  feed  character 


command  buffer 

read  data  buffer 

write  data  buffer 

board  name  buffer 

board  or  device  name  buffer 

file  name  buffer 


CLS 

DIM  W( 1500), RRR( 1500) 

DIM  RDSS$(1024) 

SCREEN  0 

COLOR  7,1,3 

OPEN  "C: DATA. DATA"   FOR  INPUT  AS  #1 

BEEP:  CLS 
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4100 

REM  *      BEGIN  MAIN  PROGRAM             * 

4120 

REM  *****************-{(**->'<->''****->'f**'^**->'<******* 

4140 

LOCATE  5,8: PRINT"   DATA  ACQUISITION     " 

4160 

REI*1  *Vc**Vf?VVf***iV>V***iVVf*iV*^V>V*>V**iV*Vf*iV***iV*iV**********Vc***->V*****'5V'!V'>V 

4180 

REM  *         T$  =  TARGET  IDENTIFICATION 

* 

4200 

REM  *                     W$  =  WAVEFORM  TYPE 

* 

4220 

REM  *         D$  =  DATE 

* 

4240 

REM  *                     N   =  NUMBER  OF  SUB -AVERAGE  /  WAVEFORM  OR  COUNT 

Vc 

4260 

REM  *         T2  =  DPO  SAMPLING  INTERVAL  (nsec) 

* 

4280 

REM  *         SI  =  DPO  MAXIMUM  DPO  VERTICAL  SCALF.(M) 

* 

4300 

REM  *         ND  =  NUMBER  OF  DATA  BLOCKS 

* 

4320 

REM  *         PI  =  NUMBER  OF  POINT  /  WAVEFORM 

* 

4340 

REM  *                     Fl  =  DPO  TIME  WINDOW 

* 

4360 

REM  'V        DSE$  =  FILE  NAME  OF  THE  DATA  OUTPUT  FILE 

* 

4380 

REM  *        DEV$  =  NAME  OF  THE  HPDEV  IN  THE  CONFIGURATION  FILE* 

4400 

REM  *        DSO%  =  THE  NUMBER  ASSOCIATED  WITH  THE  DEVICE 

4420 

REM  *        VEW$  =  STRING?  VARIABLE  FOR  SENDING  TO  HP54120 

* 

4440 

REM  *      DUMMY?  =  TIME  SET 

* 

4500 

REI>J  iVVr*VciV^VV«V*-)ViVVfVfiViV*-)VVfiV^VTV?V^VVf*iV*-!V^ViV***Vr*iV-iViV*Vf**A^V**iV***Vr*Vr*-)V*>V*** 

4520 

REM 

4540 

INPUT^n,Pl 

4560 

PRINT  "NUMBER  OF  POINT ", 

Fl 

4580 

INPUT/n,N 

4600 

PRINT  "NUMBER  OF  SUB -AVERAGE  (count) " 

N 

4620 

INPUT/n,ND 

4640 

PRINT  "NUMBER  OF  DATA  BLOCK " 

ND 

4660 

INPUT#1,T$ 

4680 

PRINT  "TARGET  ID " 

,T$ 

4700 

INPUT#1,D$ 

4720 

PRINT  "DATE " 

D? 

4740 

INPUTS 1,  T2 

4760 

PRINT  "DPO  SAMPLING  INTERVAL " 

,  T2 

4780 

INPUT//1,  W$ 

4800 

PRINT  "WAVEFORM  TYPE ";  W$ 

4820 

INPUT#1,F1 

4840 

PRINT  "DPO  TIMEWINDOWC  NSEC  ) ";  Fl 

4860 

INPUTy/l,  SI 

4880 

PRINT  "DPO  VERTICAL  SCATiE(M) ";  SI 

4900 

G1=S1 

4920 

INPUT^'U,  DSE$ 

4940 

PRINT  "DATA  FILE  NAME  OUTPUT ";  DSE$ 

4960 

INPUT/n,DUMMY$ 

4980 

PRINT  "AUTO  TIME  SET ";  DUMMY? 

5000 

INPUT#1,DELAY1 

5020 

PRINT  "DELAY  (  NSEC  ) ";DELAY1 

5040 

CL0SE#1 

5060 

PRINT  " " 

5080 

COLOR  7,1:  PRINT  "  PRESS  THE  FOLLOWING  KEY  " 

5100 

PRINT  "  Y    (  CHANGE  THE  DATA.  DAT  ) " 

5120 

PRINT  "  N    (NO  CHANGE  ) " 

5140 

V$=INKEY$:  IF  V$=""  THEN  5140 

5160 

IF  V$="Y"  OR  V$  ="y"  THEN  11260 

5180 

IF  V$="N"  OR  V$  ="n"  THEN  5220 

5200 

GOTO  5100 

5220 

DEV$="HPDEV" 

5240 

REM  VfsVVfV-VfV-y.-Vf^VfVfAVMVVrVfyttVVcvVAVfVfVfA^VfVryfAVf^A-'c*****- 

>c-fc 
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5260  REM  *  OPEN  DEVICE  AND  RETURN  THE  UNIT  * 
5280  REM  *  DESCRIPTOR  ASSOCIATED  WITH  THE  GIVEN  * 
5300  REM  *  NAME  * 

5320  REM  'VVnWfiV')V'>V')ViV**iV'>V***'jV**Vr*Vc******V«VVriV***iVVf**** 

5340  CALL  IBFIND(DEV$ ,DSO%) 

5360  IF  DSO%  <  0  THEN  PRINT  "ERROR  IN  IBFIND" 

5400  REM  *  CLEAR  SPECIFIED  DEVICE  * 

5420  REM  *****************i(*******>'f*********-f'***** 
5440  CALL  IBCLR(DSO%) 

5460  REM  *'V*Vc**iV*'>V'>V>V*'5V****************iV*'5V**'A''>'r**'A'*'sV 

5480  GOSUB  6080 

5500  CLS'BEEP 

5520  PRINT  W$+"  DATA  STORAGE  COMPLETE" 

5540  PRINT  "  MEASURE  OTHER  WAVEFORM  TYPE(Y/N)  " 

5560  Q$=INKEY$: IF  Q$="  "  THEN  GOTO  5560 

5580  IF  Q$="N"  OR  Q$="n"  THEN  GOTO  5640 

5600  IF  Q$="Y"  OR  Q$="y"  THEN  GOTO  4040 

5620  GOTO  5560 

5640  STOP 

5660  REM"''^''''******'''''"''''''***''' *'"''''''**'''■**'''''''■*************'''''* 

5680  REM*        SUBPROGRAM  * 

5700  REM*  * 

5740  REM  *  TEXT  REFERENCE  * 

5760  REM  *  "GURU  II  GPIB  USER'S  RESOURCE"     * 

5780  REM  *  "UTILITY  FOR  IBM  PC"  * 

5800  REM  *  "TEKTRONIX  1986"  * 

5820  REM  'V7ViVVr*iV'5ViVVfVfVri'rVriVy-iV'5ViVVf7ViV7ViViVt^iVT^iVVfVr'sV'iVi'r*'jV'5V'A'')VVc*')V 
5840  REM  ^'ViV^V^VVriViVVriVVrVc^ViV^VVr-.ViViVVriVVc-iVVf'jViV'jViVVrVf'A-Vf-'-'sV^iV^^'iV'jViV'jV 

5860  REM  *  WRITE  DATA  FROM  STRING  SUBROUTINE  * 
5880  REM  'V-v*-v?vv«v-V'>v-vvc-.V7V*-V')V****iV**7V**Vf7V*Vc****TV**')Sr*7ViV 

5900  CALL  IBWRT(DSO%,VEW$) 
5920  RETURN 

5940  REM  'VVr'5ViV'>VVfiV'>VV-Vf7V7V'>V'5V'5'riViV'!ViV7V7VVr7ViVyciVVfypiV')ViV*5lriV*iV**iV5V'sV 

5960  REM  *  TEXT  REFERENCE  * 

5980  REM  *  "GETTING  STARTED  GUIDE"  * 

6000  REM  *  "HP  54120T  DIGITIZING  "  * 

6020  REM  *  "OSCILLOSCOPE"  * 

6040  REM  *  "HP  COMPANY  1987"  * 

6060  REM  'V**'>V*Vr*Vr*'>ViV'>V>V'>ViViV'>VVf V-Vf ■>V'jV*yf ******'>V'5V-)V'A-*'>V'5V->V*'sV'!V 

6080  REM  "*********"***"""*"*""*■**•«■***************■ 

6100  REM  *        SETUP  THE  DISPLAY  * 

6120  REM  ***********yfyc*VcyfV-Vc*-vy-VrVr*Vfyfyf ************** 

6140  VEW$="*RST" 

6160  GOSUB  5840 

6180  REM  ***************************************** 

6200  REM  *  CRT  FORMAT  1  PROVIDE  ONE  DISPLAY  AREA  * 

6220  REM  *  AND  USE  8  DIVISIONS  FOR  THE  FULL      * 

6240  REM  *  SCALE  RANGE  * 

6260  REM  ***************************************** 

6280  VEW$=":  DISPLAY:  GRATICULE  FRAME; FORMAT  l" 

6300  GOSUB  5840 

6320  REM  " ****-'f ******'^r**************************** 

6340  REM  *  TURN  TDR  STEP  GENERATOR  ON  * 
6360  REM  *************Vf ****************Vryf *VcVf *****vr 
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6380  VEV$=": NETWORK:  REFLECTION:  PRESET" 

6400  GOSUB  5840 

6440  REM  *   PAGE  6-13  TURN  ON  CHANNEL  4         * 

6480  VEW$=":VIEW  CHANNEL4" 

6500  GOSUB  5840 

6520  REM  'V**vt-V'}VVr***?v****Vf*y«v**iViV***iViViViV**'>v***Vf**** 

6540  REM  *    PAGE  6-4  TURN  OFF  CHANNELl         * 

6560  REM  *Vf*VfV-**'jV'>VVfVr'jVyr'>VVfVfiV'>Vyf*A'A'>VVriV*?VVf:V'>V*?V'/r**Vr:'fyrTV'A'* 

6580  VEW$=": BLANK  CHANNELl" 

6600  GOSUB  5840 

6640  REM  *   SET  UP  CHANNEL  RANGE  AND  OFFSET     * 

6660  REM  *   PAGE  17-6  * 

6680  REM  *   TO  BE   20  NANOSECONDS  * 

6720  VEW$=":TIMEBASE:  RANGE  20  NS" 

6740  GOSUB  5840 

6760  REM  PRINT  ''''f^'ViViV^v^iv^VivVfiVVr^ViViVVf^Aivvr-sVTViV******^?***************^^***'' 

6780  REM  PRINT  "SET  THE  VERTICAL  SCALE  RANGE  FROM  8  TO  640  MV     " 

6840  VERT=S1 

6860  VEW$=":CHANNEL4: RANGE  "+STR$(VERT)+"m" 

6880  GOSUB  5840 

7040  REM  *       TURN  ON  THE  GRATITUDE  * 

7080  VEW$=": DISPLAY:  GRATICULE  GRID" 

7100  GOSUB  5840 

7120  VEW$=":  ACQ: BAND  LOW" 

7140  GOSUB  5840 

7180  REM  *       AVERAGING  ON  * 

7220  VEW$=": ACQUIRE: TYPE  AVERAGE" 

7240  GOSUB  5840 

7280  REM  *    ASKING  FOR  BANDWIDTH  * 

7300  REM    'V**'>V';VV--VV-':V'>ViVVr^ViVVr?ViVvVVr'5VVr7V->'5V*7ViVVf ******'A''>V'>V>V'5ViV'5V 

7320  PR  I  NT    ' '  i"- "  ^fVf  Vj'^rV-VcV-V.-Vr  V-VcVf  VrVr-VV-^Vj  Vr  V--V-V-V*yr-!'ryf  ■A'^^iVVcVcVf  Vf ' ' 

7340  PRINT  "THE  DEFAULT  BANDWIDTH  IS  12. 4GHZ     " 

7360  PRINT  "WOULD  YOU  LIKE  TO  CHANGE  TO  20  GHZ  (Y/N)" 

7380  Q$=INKEY$: IF  Q$="  "  THEN  GOTO  7380 

7400  IF  Q$="N"  OR  Q$="n"  THEN  GOTO  7500 

7420  IF  Q$="Y"  OR  Q$="y"  THEN  GOTO  7460 

7440  GOTO  7380 

7460  VEW$=": ACQ: BAND  HIGH" 

7480  GOSUB  5840 

7500  RD$=SPACE$(255) 

7520  VEW$="ACQUIRE:  COUNT  "+STR$(N) 

7540  GOSUB  5840 

7560  VEW$=": ACQUIRE: POINTS  "+STR$(P1) 

7580  GOSUB  5840 

7600  REM  *Vfyf*VfVrVf?VVfV-yf*Vc*yfVryfVfyfVfVryc'Vf'^Vf***yrVr**Vfyc-'A'*':V**'jV 

7620  REM  *         SET  TIME  DELAY  * 

7640  REM  *  * 
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7660  REM  *****iv************iv ********************* 

7680  REM 

7700  REM 

7720  REM 

7740  DELAY=DELAY1 

7760  DUM$=STR$( DELAY) 

7780  VEW$=":TIM:DEL  ^'+DUM$+"NS" 

7800  GOSUB  5840 

7820  REM 

7840  REM 

7860  REM 

7880  REM 

7900  REM 

7920  REM 

7940  REM 

7960  REM  **************************************** 

7980  REM  HP  54120T  SENT  THE  LAST  BYTE  RESPONSE  "LF"  ASH  DECIMAL  10 

8000  VEW$=":  SYSTEM:  HEADER  OFF; : EOI  ON" 

8020  GOSUB  5840 

8040  PR I NT  ' '  **************************************************" 

8320  REM  **************************************** 

8340  PRINT  "I  AM  GOING  TO  DIGITIZE  CHANNEL(Y/N)" 

8360  Q$=INKEY$:  IF  Q$="  "  THEN  GOTO  8360 

8380  IF  Q$="Y"  OR  Q$="y"  THEN  GOTO  8420 

8400  GOTO  8360 

8420  PRINT  TIME$ 

8560  REM  ***************************************** 

8600  REM  *    GET  THE  NUMBER  AND  STORE  IN  THE    * 

8620  REM  *    FILE  * 

8640  REM  ***************************************** 

8660  DSEW$="C:C" 

8680  FOR  1=1  TO  ND 

8700  IF  I<10  THEN  DSFILE$=DSEW$+RIGHT$( STR$( I) , 1)+".  WFM" 

8720  IF  I>9  THEN  DSFILE$=DSEW$+RIGHT$(STR$( I) ,2)+".  WFM" 

8740  REM 

8760  REM  ***************************************** 

8780  REM  *  DIGITIZE   CHANNEL  4  * 

8800  REM  ***************************************** 

8820  VEW$=": DIGITIZE  CHANNEL4" 

8840  REM 

8860  GOSUB  5840 

8880  VEW$=": WAVEFORM:  SOURCE  WMEM0RY4;  FORMAT  ASCII" 

8900  GOSUB  5840 

8920  VEW$=": WAVEFORM:  DATA?" 

8940  GOSUB  5840 

8960  REM 

8980  CALL  IBRDF(DSO%,DSFILE$) 

9000  CLOSE/n 

9020  NEXT  I 

9040  REM  ***************************************** 

9060  REM  *   FINISH  DIGITIZING  * 

9080  REM  ****************VcVf*********************Vf* 

9100  PRINT  TIME? 
9120  DSFILE$="HEADER" 
9140  VEW$=":WAV:PRE?" 
9160  GOSUB  5840 
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9180  CALL  IBRDF(DSO%,DSFILE$) 

9200  PRINT  "HEADER  COMPLETE" 

9220  GOSUB  10200 

9240  PRINT  "'V'}V-v*Vf>v*Vc'>vvfiVVf!V**iV?v*'>v>v*" 

9260  PRINT  "WOULD  YOU  LIKE  TO  GET  A  HARD  COPY?(Y/N)" 

9280  Q$=INKEY$:  IF  Q$="  "  THEN  9280 

9300  IF  Q$="N"  OR  Q$="n"  THEN  GOTO  10080 

9320  IF  Q$="Y"  OR  Q$="y"  THEN  GOTO  9360 

9340  GOTO  9280 

9360  REM  'VVoViVVrVoV^^yr^V^V^VVciViViV^^ViV^V^V^ViV^V^Vyr^^A^V^Vc******* 

9380  REM  *   PLOT  THE  LAST  WAVEFORM  * 

9400    RE  M    'ViV^^VriVVf VrViriVVf Vf ■jVVc'jV'jVVcVciVVc'jVVc'Vf VriVVr A''A'iV;lrVr^VVc*ilr*^V'^')<rVoV 

9420  VEW$=": HARD:  SOUR  FACT,WMEM4" 

9440  GOSUB  5840 

9460  PRINT  "PLOT  ?  (Y/NV' 

9480  Q$=INKEYS:  IF  Q$="  ^'  THEN  9480 

9500  IF  Q$="N"  OR  Q$="n"  THEN  GOTO  9600 

9520  IF  Q$="Y"  OR  Q$="y"  THEN  GOTO  9560 

9540  GOTO  9480 

9560  VEW$=":PLOT?;" 

9580  GOSUB  5840 

9570  GOTO  9780 

9600  REM  *■?'{■  ?VVriV-;VV.-iV>VVc*')ViV*Vf iV^ViV^-iWc VnVVr** *VrVr***yfVc***TiV?V** 

9620  REM  *   PRINTING  * 

9660  VEW$=":HARDCOPY:  PRINTER  DEFAULT" 

9680  GOSUB  5840 

9700  VEW$=":HARDCOPY:PAGE  AUTOMATIC" 

9720  GOSUB  5840 

9740  VEW$=":  PRINT?;" 

9760  GOSUB  5840 

9780  B0AD$="GPIB0" 

9800  CALL  IBFIND(B0AD$,B0A%) 

9820  CALL  IBSIC(BOA%) 

9840  VEW$="?Q!" 

9860  CALL  IBCMD  (BOA%,VEW$) 

9880  V%=1 

9900  CALL  IBGTS  (BOA%,V%) 

9920    PRINT    "'V'jWfiV^iVyrVciVVc^VVv-.VVciViVVciViV^iVTVVfiVVf^V^A^VfiViViV^^iViViV" 

9940  PRINT  "PLEASE  WAIT  UNTIL  THE  PLOTTER  HAD  FINISHED" 

9960  PRINT  "PLEASE  RELEASE  THE  DPO  TO  LOCAL(Y)" 

9980  Q$=INKEY$: IF  Q$="  "  THEN  9280 

10000  IF  Q$="Y"  OR  Q$="y"  THEN  GOTO  10040 

10020  GOTO  9980 

10040  CALL  IBLOC(DSO%) 

10060  CALL  IBLOC(BOA%) 

10080  RETURN 

10100  REM  *****''"f*'5'f****'/f***«'f*'/f*'A'*Vf****VcyrVfVf ********* 

10120  REM  SUBROUTINE  READ  DATA  TO  STRING  REF.  GURU  D-47 

10140  CALL  IBRD(DSO%,RD$) 

10160  RETURN 

10180  REM  *'**"A"'^<'**''«"A'*'^f*'/r''''^<'*****Vc**********Vf ******** 

10200  REM  *  PROGRAM  TO  ARRANGE  FORMAT  FILE   * 

10220  REM  *  SCALE  THE  DATA                  * 

10240  REM  *  * 

10260  REM  *  * 
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10280 
10300 
10320 
10340 
10360 
10380 
10400 
10420 
10440 
10460 
10480 
10500 
10520 
10540 
10560 
10580 
10600 
10620 
10640 
10660 
10680 
10700 
10720 
10740 
10760 
10780 
10800 
10820 
10840 
10860 
10880 
10885 
10890 
10892 
10896 
10897 
10898 
10899 
10900 
10901 
10902 
10903 
10904 
10910 
10920 
10940 
10960 
10980 
11000 
11020 
11040 
11060 
11080 
11100 
11120 
11140 


REM  ***'A">'"V'jV*'jV'>v 
OPEN  "C: HEADER 
INPUT#1,A%,B%, 


"FORMAT 

"TYPE 

"POINT 

"COUNT 

"XINCRE 

"XOR 

"XREF 

"YINC 

"YOR 

"YREF 

"RANGE 


FOR  INPUT  kSin 
C% , D % jE,F,G/4,H,  I  ,J,K 
■"  A% 


B% 

C% 

D% 

E 

F 

G% 

H 

lY 

J 

K 


PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

CLOSE/n 

REM   DSE$="C:  C" 

FOR  K=l  TO  ND 

IF  K<10  THEN  DSFILE$=DSEW$+RIGHT$(STR$(K),l)+".wfm" 

IF  K>9  THEN  DSFILE$=DSEW$+RIGHT$(STR$(K) ,2)+".  wfm" 

IF  K>99  THEN  DSFILE$=DSEW$+RIGHT$(STR$(K) ,3)+".  wfm" 

OPEN  DSFILE?  FOR  INPUT  AS  in 

REM 

REM 

REM 

REM 

REM 

REM 

FOR  1=1  TO  C% 
INPUT  #1,WFM 
W(I)-=-  ((WFM-J)*H+IY) 
NEXT  I 
SUMM=0.  0 

FOR  1=1  TO  C% 
SUMM=SUMM+W(I) 

NEXT  I 

REM  *  REMOVE  MEAN  * 

SMEAN=SUMM/C% 
PRINT  "MEAN  ="; SMEAN 

FOR  1=1  TO  C% 

RRR(I)=W(I)-SMEAN 

NEXT   I 
CLOSE/n 

DSFILE$=DSE$+".WFM" 
OPEN  "A",/n,DSFILE$ 
PRINTS 1,T$ 
PRINT/n,W$ 
PRINTS 1,D$ 
PRINTS 1,DUMMY$ 
PRINT/n.Dro 

PRINT?'n,USING  "iHHh  W;  Fl 
PRINTyn.ND 
PRINTS 1, PI 
FOR   1=1  TO  C% 
PRINT/n,USING"#//.  imU  ";  RRR(I) 
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11160 
11180 
11200 
11220 
11221 
11240 
11260 
11280 
11300 
11320 
11340 
11360 
11380 
11400 
11420 
11440 
11460 
11480 
11500 
11520 
11540 
11560 
11580 
11600 
11620 
11640 
11660 
11680 
11700 
11720 
11740 
11760 
11780 
11800 
11820 
11840 
11860 
11880 
11900 
11920 
11940 
11960 
11980 
12000 
12020 
12040 
12060 
12080 
12100 
12120 
12140 
12160 
12180 
12200 
12220 
12240 


NEXT  I 
REM 

CLOSED 1 
NEXT  K 
BEEP: BEEP 
RETURN 


REM  * 
REM  * 
REM  * 
REM  * 
REM  * 
REM  * 
REM  * 
REM  * 
REM  * 
REM  * 
REM  * 


PI 

N 

ND 

T$ 

D$ 

T2 

W$ 

Fl 

SI 


it 

* 
* 
* 
* 
* 
* 

* 

TS$=SPACE$(64) 

DUMMY$=SPACE$(64) 
KEY  OFF:  SCREEN  0,0,0 
OPEN  "DATA. DAT"  FOR  INPUT  AS  in 

INPUT//1,P1:  INPUT//1,N:  INPUT#1,ND:  INPUTyn,T$:  INPUT#1,D$:  INPUT^/l 
INPUT#1,W$:  INPUT//1,F1:  INPUT#1,S1:  INPUTS  1  ,DSE$:  INPUT#1,DUMMY$ 
INPUTS/ 1,  DELAY  1:CL0SE#1    .__^ 
COLOR  7,1,3:  CLS  ^~^ 

X=4:  Y=16 
PRINT"PROGRAM  DATA  ACQUISTION 


=  NUMBER  OF  POINTS 
=  NUMBER  OF  SUBAVERAGES 
=  NUMBER  OF  DATA  BLOCKS 
=  TARGET  IDENTIFICATION 
=  DATE  OF  MEASUREMENT 
=  DPO  SAMPLING  INTERVAL 
=  WAVEFORM  TYPE 
=  DPO  TIME  WINDOW 
=  DPO  MAXIMUM  SCALE 

DSE$  =  OUTPUT  DATA  FILE  NAME 

DUMMY $  =  DUMMY  $ 


,T2 


LOCATE  X,Y+15:  PRINT 
LOCATE  X+2,Y: PRINT 
LOCATE  X+3,Y: PRINT 
LOCATE  X+4,Y:  PRINT 
LOCATE  X+5,Y: PRINT 
LOCATE  X+6,Y: PRINT 
LOCATE  X+7,Y: PRINT 
LOCATE  X+8,Y:  PRINT 
LOCATE  X+9,Y:  PRINT 
LOCATE  X+10,Y: PRINT 
LOCATE  X+11,Y: PRINT 
LOCATE  X+12,Y:  PRINT 
LOCATE  X+13,Y: PRINT 
LOCATE  X+14,Y:PRINT 
Q$=INKEY$: IF  Q$="  " 


PRESS  THE  FOLLOWING  KEY" 
<  F  >  FILE  NAME  FOR  DATA  OUTPUT. . . ! 
NUMBER  OF  POINTS 


P 
A 
B 
T 
D 
W 
K 


>  NUMBER  OF  SUB  AVERAGES. 

>  NUMBER  OF  DATA  BLOCKS. . 

>  TARGET  ID.  ($) 

>  DATE   ( $ ) 

>  WAVEFROM  TYPE 

>  DPO  TIME  WINDOW  

S  >  DPO  SAMPLING  INTERVAL. 

R  >  RUN  THE  PROGRAM 

U  >  DPO  MAXIMUM  VERT  MV. . . 

X  >  AUTO  TIME  SET 

G  >  DELAY  TIME 


DSE$ 

PI 

N 

ND 

T$ 

D$ 

W$ 

Fl 

T2 

SI 

DUMMY $ 
38 


IF  Q$="G' 
IF  Q$="F' 
IF  Q$="P' 
IF  Q$="A' 
IF  Q$="B' 
IF  Q$="T' 
IF  Q$="D' 
IF  Q$="W' 
IF  Q$="S' 
IF  Q$="K' 
IF  Q$="R' 
IF  Q$="U' 


OR  Q$= 
OR  Q$= 
OR  Q$= 
OR  Q$= 
OR  Q$= 
OR  Q$= 
OR  Q$= 
OR  Q$= 
OR  Q$= 
OR  Q$= 
OR  Q$= 


'k" 
II  It 

r 

OR  Q$="u" 


•f" 

t  II 

P 

•a" 


THEN  12000 
THEN  12900 
THEN  12300 
THEN  12340 


THEN  12560 

'b"  THEN  12600 

't"  THEN  12640 

'd"  THEN  12680 
t  II 
w 

I   IT 

s 


THEN  12700 
THEN  12740 
THEN  12780 
THEN  4040 
THEN  12820 
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12260 
12280 
12300 
12320 
12340 
12360 
12380 
12400 
12420 
12440 
12460 
12480 
12500 
12520 
12540 
12560 
12580 
12600 
12620 
12640 
12660 
12680 
12700 
12720 
12740 
12760 
12780 
12800 
12820 
12840 
12860 
12880 
12900 
12920 
12940 
12960 
12980 
13000 
13020 
13040 
13060 


IF  Q$="X"  OR  Q$="x"  THEN  12860 

GOTO  11720 

PRINT" INPUT  FILE  NAME($)  EX. . C: F" 

INPUT  DSE$:GOTO  12940 

CLS:  LOCATE  4,8: PRINT"NUMBER  OF  POINT 

PRINT  "A).  .  .  1024 POINT. 


512 POINT. 

256 POINT. 

128 POINT. 


Q$="a" 
Q$="b" 


THEN 
THEN 
THEN 
THEN 


Pl=1024 
Pl=512 
Pl=256 
Pl=128 


PRINT  "B). 

PRINT  "C). 

PRINT  "D). 

Q$=INKEY$:  IF  Q$=""  THEN  12440 

IF  Q$="A"  OR 

IF  Q$="B"  OR 

IF  Q$="C"  OR 

IF  Q$="D"  OR  Q$="d" 

GOTO  12940 

PRINT" INPUT  NUMBER  OF  SUB -AVERAGE" 

INPUT  N:GOTO  12940 

PRINT  "INPUT  NUMBER  OF  DATA  BLOCK" 

INPUT  ND  :GOTO  12940 

PRINT" INPUT  TARGET  ID($)" 

INPUT  T$:GOTO  12940 

PRINT"SET  DATE": D$=DATE$: GOTO  12940 

PRINT" INPUT  WAVEFROM  TYPE($V' 

INPUT  W$: GOTO  12940 

PRINT" INPUT  DPO  SAMPLING  INTERVAL" 

INPUT  T2:G0T0  12940 

PRINT"DPO  TIME  WINDOW         " 

INPUT  Fl:  GOTO  12940 

PRINT" INPUT  DPO  VERTICAL  SCALE" 

INPUT  SI: GOTO  12940 

PRINT" AUTO  TIME" 

DUMMY$=TIME$:  GOTO  12940 

PRINT" INPUT  TIME  DELAY" 

INPUT  DELAY1:G0T0  12940 

OPEN  "DATA. DAT"  FOR  OUTPUT  AS#1 

PRINT#1 ,  PI:  PRINT#1  ,N:  PRINT#1  ,ND:  PRINT#1  ,T$:  PRINTj^/l  ,D$ 

PRINT#1,T2:  PRINT#1,W$:  PRINT/n  , USING  "###.M";F1:  PRINT#1,S1 

PRINTS 1,DSE$: PRINT#1 , DUMMY?: PRINT#1, DELAY 1 

CL0SE#1  :  CLS 

REM 

GOTO  11700 
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APPENDIX  C.      CHARACTERISTICS  OF  LOW  LOSS  CABLE 


Table  21.      CHARACTERISTICS  OF  PORT  RETURN  CABLES 


Impedance 

50  ohms 

Capacitance 

26  pF/ft 

Time  delay 

1.2ns;ft 

Velocity  of  propagation 

85%  of  light  velocity 

Jacket  withstand 

1.0  kV 

Center  conductor 

15  AWG  sohd 

Minimum  bend  radius 

1  inch 

Dielectric  constant 

1.4 

Dielectric  withstand 

I.OkV 

Table  22.      INSERTION  LOSS 


0-4  GHz 

less  than  0.36  dB 

4-8  GHz 

less  than  0.5   dB 

8-12  GHz 

less  than  0.6   dB 

12-16  GHz 

less  than  0.72  dB 

16-18  GHz 

less  than  0.8   dB 

1S-26.5GHZ 

less  than  1.0  dB 
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